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An  error  was  found  in  the  application  of  the  formula  providing  the  L 
and  L1  adjusting  terms  in  the  second  approximation  stage,  as  carried 
out  in  the  previously  issued  report.  This  mistake  affects  most  of 
the  subsequent  results  to  a small  degree,  but  in  no  case  are  the 
actual  duct  coordinates  changed  by  as  much  as  one  percent.  The  cor- 
rections have  been  meticulously  made,  however,  beginning  on  page  48, 
in  order  not  to  detract  from  the  educational  value  of  the  report,  but 
nothing  of  any  consequence  in  regard  to  design  can  be  attributed  to 
the  slight  numerical  changes  which  have  resulted  from  these 
recomputations. 
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SUMMARY 


Several  numerical  examples  have  been  worked  out  in  detail,  by  applica- 
tion of  a theoretical  procedure  developed  in  Part  I of  this  study,  in  order  to 
determine,  consistent  with  the  accuracy  inherent  in  the  use  of  linear  theory, 
what  the  optimal  contour  of  an  annular  duct  must  be  if  it  is  to  produce  min- 
imum external  wave  drag  in  supersonic  flow,  under  the  stipulations  that  either 
(1)  the  area  enclosed  between  the  sought  contour  and  its  inner  cone-like  wall 
is  to  have  a set  constant  value,  measured  in  any  meridional  plane  through  the 
duct's  axis,  or  (2)  the  volume  enclosed  between  the  surface  of  revolution  swept 
out  when  the  sought  meridional  contour  is  rotated  about  its  axis  and  the 
frustum-of-a-cone  inner  surface  is  to  be  a constant  given  value. 

The  specific  values  that  have  been  selected  for  the  parameters  which 
describe  the  general  features  of  the  free-stream  flow  and  the  duct's  overall 
geometry  are  such  as  to  ensure  enough  generality  so  as  to  cover  practically 
all  cases  of  interest  in  practical  applications.  In  all  the  cases  examined  it 
is  found  that  the  coordinates  of  the  best  ducts,  which  correspond  to  the 
various  perimetral  conditions  that  are  involved,  turn  out  to  be  scarcely  dis- 
tinguishable from  the  contour  which  would  be  obtained  in  each  specific  case 
by  making  use  of  the  simple  assumption  that  the  pressures  created  on  the  duct 
are  directly  related  to  the  local  slope  of  the  planar  elements  according  to 
the  dictates  of  Ackeret's  two-dimensionally  valid  formula.  This  far-reaching 
conclusion  greatly  simplifies  the  problem  of  the  designer  in  determining  such 
optimum  duct  contours  in  practice,  but  one  must  not-  lose  sight  of  the  fact 
that  the  local  pressures  obtained  by  use  of  the  Ackeret  formula  are  not  go- 
ing to  be  just  minor  variants  from  the  actual  pressures  which  would  be 
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obtained  by  use  of  the  truly  three-dimensional  approach  described  in  the  theory 
now  under  examination.  The  amounts  that  the  local  pressures  on  the  various 
ducts  differ,  between  the  two-  and  three-dimensional  treatments,  are  not 
worked  out  in  all  detail,  but  a comparison  of  the  general  harmonic,  terms, 

which  go  to  make  up  the  total  axial  components  of  the  induced  velocities  at 

/ 

the  duct's  surface,  shows  that  the  terms  of  small  order  are  significantly 
different  in  the  two  contrasting  approaches. 
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j I.  INTRODUCTION 

(1) 

, The  theory  developed  In  Part  I of  this  study  is  now  to  be  worked  out  in 

numerical  detail  for  the  general  case  defined  by  the  values 

pB  = 0.25  and  £,=  Z (l) 

where  B is  the  cotangent  of  the  Mach  angle  and  p is  the  semi-vertex  angle  of  the 
frustum  of  a cone  which  constitutes  the  inner  surface  of  the  cone-like  central  duct 
of  the  body  of  revolution  being  sought-  In  addition,  represents  a 

"reduced"  length  for  the  duct  inasmuch  as  it  stands  for  the  axial  length  of  thedueti 
measured  in  abscissa  units  which  may  be  defined  as  the  maximal  distance  out  ahead 
of  the  duct,  and  on  its  axis,  where  a small  disturbance  must  be  located  in  order 
that  it  should  be  felt  only  downstream  of  the  position  of  the  duct's  lip  at  radial 
distances  equal  to  or  larger  than  the  one  defining  the  size  of  the  duct's  mouth* 

The  external  contour  of  the  body  of  revolution  having  such  a central  duct 
is  then  to  be  determined  in  such  wise  as  to  give  minimum  drag  tinder  one  or  more  of 
several  perimetral  conditions  that  are  further  imposed  on  the  geometry  of  the 
sought  missile  shape.  The  further  conditions  which  are  thought  to  be  of  some  prac- 
tical consideration  are  of  the  following  types: 

(a)  No  perimetral  condition  imposed,  except  to  assume  that  the  external  contour 

actually  closes  without  a blunt  base;  i*e.,  it  is  merely  assumed  that  the  "end- 
points" N-^  and  $2  at  the  lip  and  at  the  tail  end  of  the  sought  duct  contour, 

I 

respectively,  are  to  belong  to  the  frustum-of-a-cone  section  which  constitutes 

the  inner  duct  whose  large  end  and  small  end  have  radii  in  the  ratio 

— = 1 + p * ^ ; where,  of  course,  Ri  is  now  equal  to  the  radial  coordinate 
K1  R1 

of  the  N]_  lip-point  and  R2  is  equal  to  the  radial  coordinate  of  the  N2  tail- 
point  of  the  external  contour-  This  example  is  worked  out  for  two  specific 
free-stream  flow  conditions,  defined  as 
Normalizing  Mach  No.  Case: 

M^  = i~2,  so  B = / Mj'  -1  - 1,  and  thus  p = 0-25  and  L/R-^  = 3 
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Small  Conical  Diffusor  Angle  Case; 


Mco=  1 9436,  so  B = / mJ -1  = 1.6667,  and  thus  0 = 0.15  and  b/\  = 5. 

(This  type  of  duct  might  better  be  designated  merely  the  High  Fineness  Ratio 
Case. since  the  diffusion  process  taking  place  in  the  internal  flow  is  entirely 
divorced  in  this  study  from  the  flow  over  the  external  contours  being  sought; 
nevertheless,  there  is  a certain  amount  of  sanction  to  be  found  for  this  usage 
in  the  parlance  of  jet  engine  designers  who  are  most  interested  in  the  avail- 
able SDace  in  the  nose  even  though  the  customary  use  of  a spike  central  body 
will  normally  decrease  the  cross-sectional  area  for  the  internal  flow  imme- 
diately behind  the  entrance). 

(b)  In  addition  to  the  closure  condition,  it  is  postulated  that  the  area  enclosed 

between  the  sought  contour  and  the  x-axis  is  to  be  equal  to  the  area  which 

would  be  enclosed  between  this  x-axis  and  the  straight  line  segment  joining 
points  and  N2  (this  case  will  be  designated  as  '"No  Increase  in  Area  Over 
that  Occupied  by  Frustum-of-Cone  Basic  Shape”) > Of  course,  the  inner  duct 
will  have  to  have  some  other  shape  than  the  frustum  of  a cone  delimited  by 
the  Ni  and  N2  end  points,  but  the  flow  through  the  inner  duct  is  not  being 
considered  in  any  of  these  studies.  Same  free-stream  flow  conditions  as  in  (a). 

(c)  In  addition  to  the  closure  condition,  it  is  postulated  that  the  area  enclosed 

between  the  sought  contour  and  the  chord  line  joining  the  Nj.  and  N2  endpoints 

is  to  be  equal  to  the  area  which  would  be  enclosed  between  this  chord  line 
between  Nj_  and  N2  and  lying  under  the  curve  with  equation 

where  the  system  of  Cartesian  coordinates  employed  here  in  a meridional  plane 
through  the  body  axis  is  so  orientated  that  the  x-axis  and  body-axis 
coincide,  while  the  R-axis  is  taken  to  pass  through  the  lip  end-point,  Nj. 

It  may  be  noted  that  this  parabolic  contour,  used  for  sake  of  comparison,, 
is  such  that  at  the  mid-point  along  the  duct  axis  where  x/L  =0.5  the 
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"thickness"  of  the  conical  duct  fairing  is  5 % of  the  duct  length*  i.e.,  at 

R-(Rl  + 0.5  p L) 

x/L  = 0.5  it  is  found  that  AR/L j- c — - 0.05. 

The  free-stream  flow  conditions  sre  the  same  as  in  (a). 

(d)  In  addition  tc  the  closure  condition,  it  is  postulated  that  the  volume 
enclosed  between  the  solid  of  revolution  swept  out  by  the  scught  contour 
when  revolved  about  the  duct  avis  and  the  frustum  of  e cone  which  has  the 
chordline  Ni  No  as  a generatrix  shall  be  equal  to  the  volume  enclosed  between 
this  same  inner  truncated-cone  type  of  body  of  revolution  and  the  one  which 
is  swept  out  as  a result  of  a complete  rotation,  about  iie  duct's  central 
axis,  of  the  parabolic  surface  defined  according  to  the  equation 


in  the  case  of  the  Normalizing  Mach  No.:M^  = ~/~2,  or  B - 1,  and  p - 0.25, 

while  the  generating  curve  is  similarly  taken  to  be 

in  the  case  of  the  Small  Conical  Diffusor  Angle,  when  = 1.9436,  B = 1.6667, 

and  p = 0.15 

In  addition,  the  parabolic  reference  meridonal  line  defined  in  paragraph 
(c)  has  also  been  used  in  applying  this  volume-type  of  perimetral  condi- 
tion for  both  the  Normalizing  Mach  No.  Case  and  the  Small  Conical  Diffusor 
Angle  Case. 

The  procedure  to  be  followed  in  the  numerical  working  out  of  these 
examples  is  the  one  discussed  in  Part  I of  the  theoretical  section  of  this 
study.  Because  of  the  nature  of  the  results  found  after  making  the  numerical 
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a opli cations  for  the  examples  under  consideration  according  to  this  first 
procedure,  it  is  evident  that  there  would  be  very  little  point  in  going  on 
to  rederive  the  optimum  contours  by  use  of  the  second  of  the  two  alternative 
theoretical  approaches  presented  in  Part  II  of  the  previous  work.  As  a 
matter  of  fact,  one  may  maKe  the  following  sweeping  generalization  from 
aporaisal  ,of  the  numerical  results  calculated  in  what  follows:  as  far  as 
the  determination  of  the  optimal  external  duct  contour  is  concerned,  it  will 
require  close  observation  to  detect  any  difference  whatsoever  between  the 
best  duct  shaoes  which  result  from  use,  on  the  one  hand,  of  Eq„  (47)  of 
Part  I for  definition  of  the  relationshio  which  holds  between  the  geomet- 
rical features  of  the  sought  contour  and  the  local  pressure  coefficient, 
p— o. 

L22 , and,  on  the  other  hand,  from  use  of  the  formula  which  is  suoposed 

to  hold  strictly  merely  for  the  case  of  two-dimensional  flow  (Ackeret's  air- 
foil theory) » 

Although  this  result  does  constitute  a notable  simplification  to  the 
practical  solution  of  the  originally  set  problem,  it  may  appear  to  be  dis- 
appointing and  of  little  consequence,  especially  after  going  through  all  the 
details  required  in  the  three-dimensional  treatment.  Nevertheless,  it  should 
not  be  interpreted  as  being  of  only  minor  imoortance,  because  one  will  also  be 
able  to  observe,  on  the  basis  of  these  same  numerical  results,  to  be  presented 
below,  that  the  local  pressures  exhibit  very  sizeable  differences  between  the 
case  when  formula  (47)  of  Part  I is  employed  to  relate  the  geometrical  oroper- 
ties  of  the  duct's  contour  to  the  pressure  coefficients,  in  contrast  to  the 
case  in  which  the  Ackeret  formula  is  used. 

It  is  worth  pointing  out  in  this  regard  that  the  numerical  values  of 
PE  and  of  selected  for  purposes  of  illustration  herein,  are  such  as 
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to  make  it  probable  that  the  above-mentioned  overall  conclusion  will  still 
hold  true  for  all  other  cases  which  may  arise  in  actual  practice,  inasmuch  as 
any  cases  of  practical  interest  which  are  likely  to  be  met  with  will,  in 
general,  be  found  to  be  characterized  by  values  of  pB  and  of  Jt,  that  will 
be  smaller  in  magnitude  than  the  ones  premised  in  Eq.  (l) . 

2.  Evaluation  of  the  P Function.  Appearing  in  Eq.  (581)  of  Part  X 
The  first  step  in  carrying  out  the  computations  for  the  determination  of 
the  optimal  duct  contour  is  the  evaluation  of  the  resolvent  kernel  j1  of  the 
governing  integral  equation,  given  as  Eq.  (37)  of  Part  I.  The  expression  for 
the  resolvent  kernel  to  be  used  is 

M 


r*x)  = m - J~pg-  H%+  H (X)  -fcssJK  H <*■>  + JiTpe?  ^ ' 


(2) 


dK 


where  H(x)  = ^ , and  where  the  expression  for  K(x)  is  given  as  Eq.  (34)  of 
Part  I.  In  the  specific  case  under  study,  one  finds,  by  insertion  of  the 
numerical  value  of  pB  selected  in  Eq.  (l) , that 


HM  = d 450/6 


(3) 


and  the  moduli  of  the  complete  elliptic  integrals  El  and  appearing 

herein  are  reduced  to  just  ^ *~x  in  both  cases. 

/ 

The  values  for  H(x)  and  for  the  successively  higher  orders  of  the  iter- 
ative kernels  , H^,  , and  , as  well  as  the  pertinent  values  for 

the  resolvent  kernel  p*  itself,  have  been  determined  for  a goodly  number  of 
points  within  the  range  of  x-values  running  from  0 to  4.2,  and  the  results 
have  been  entered  in  the  appended  Table  I. 
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TABLE  I 


RESOLVENT  KERNEL  FOR  EC- (57)  OF  PART  I.  TOGETHER  WITH  ITS  INTERVENING  ITERATIVE 
KERNELS  UNDER  ASSUMPTION  THAT  BB  = 0.25 


■ 

X 

H(x) 

H(2)(x) 

*x* 

to 

H<4)(x) 

*,  v 
r w 

0 

0.38433 

0 

0 

0 

0 

0.38433 

0.1 

0.37109 

0.01409 

0.00026 

0 

0 

0.35370 

0.2 

0.35913 

0.02692 

0.00100 

0.00003 

0 

0.33583 

0.3 

0.34822 

0.03865 

0.00214 

0.00008 

0 

0.31530 

0.4 

0.33824 

0.04941 

0.00360 

0.00018 

0.00001 

0.29681 

0.5 

0.32903 

0.05932 

0.00554 

0.00033 

0.00002 

0.28001 

0.6 

0.32056 

0.06946 

0.00732 

0.00053 

0.00003 

0.26481 

0.7 

0.31266 

0.07692 

0.00948 

0.00079 

0.00005 

0.25090 

0.8 

0.30535 

0.08476 

0.01180 

0.00110 

0.00008 

0.23824 

0.9 

0.29350 

0.09205 

0.01425 

0.00148 

0.00012 

0.22659 

1 

0.29209 

0.09884 

0.01681 

0.00192 

0.00017 

0.21588 

1.2 

0.23042 

0.11563 

0.02311 

0.00309 

0.00032 

0.19348 

1.4 

0.27003 

0.13035 

0.02975 

0.00455 

0.00053 

0.17434 

1.6 

0.26071 

0.14333 

0.03660 

0.00621 

0.00080 

0.15786 

1.3 

0.25228 

0.15487 

0.04356 

0.00813 

0.00115 

0.14353 

2 

0.24462 

0.16517 

0.05058 

0.01027 

0.00158 

0.13101 

2.2 

0.23763 

0.17442 

0.05759 

0.01259 

0.00208 

0.12001 

2.4 

0.23118 

0.13276 

0.06455 

0.01509 

0.00267 

0.11026 

2.6 

0.22524 

0.19031 

0.07143 

0.01775 

0.00333 

0.10159 

2.8 

0.21974 

0.19717 

0.07822 

0.02054 

0.00408 

0.09388 

3 

0.21462 

0.20342 

0.08489 

0.02346 

0.00490 

0.08694 

3.3 

0.20755 

0.21664 

0.09805 

0.02919 

0.00655 

0.07552 

3.6 

0.20116 

0.22825 

0.11075 

0.03518 

0.00840 

0.06582 

3.9 

0.19532 

0.23350 

0.12299 

0.04137 

0.01045 

0.05748 

4.2 

0.13995 

0.24762 

0.13478 

0.04772 

0.01268 

0.05026 

p*(x)  = H(x)  - 0.89443  H^(x)  + 0.8  H^(x)  - 0.71554  H^(x) 


(5),  \ 

0.64  H (x) 
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to  make  it  probable  that  the  above-mentioned  overall  conclusion  will  still 
hold  true  for  all  other  cases  which  may  arise  in  actual  practice,  inasmuch  as 
any  cases  of  practical  interest  which  are  likely  to  be  met  with  will,  in 
general,  be  found  to  be  characterized  by  values  of  (3B  and  of  Ji  that  will 
be  smaller  in  magnitude  than  the  ones  premised  in  Eq.  (l) . 


2„  Evaluation  of  the  P Function,  Appearing  in  Eq-  (58  0 of  Part  I. 

The  first  step  in  carrying  out  the  computations  for  the  determination  of 
the  optimal  duct  contour  is  the  evaluation  of  the  resolvent  kernel  j1  of  the 
governing  integral  equation,  given  as  Eq.  (37)  of  Part  I.  The  expression  for 
the  resolvent  kernel  to  be  used  is 


r (x>  - nix)  - 


/ 


7W 


//%  + 


, w / ft)  4 
H W ~ (l-h/SBjl*-  ^ * 


H 


% • 


(2) 

dK 

where  H(x)  = ^ , and  where  the  expression  for  K(x)  is  given  as  Eq.  (34)  of 
Part  I.  In  the  specific  case  under  study,  one  finds,  by  insertion  of  the 
numerical  value  of  f3B  selected  in  Eq.  (l),  that 


H(x)  = 0.  450/6 


u+x)  E-F 

x/x-hJ.6 


(3) 


and  the  moduli  of  the  complete  elliptic  integrals  fEf  and 
herein  are  reduced  to  just  * y-  in  both  cases. 


If"  appearing 


The  values  for  H(x)  and  for  the  successively  higher  orders  of  the  iter- 
ative kernels  H^,  H^,  H^,  and  H^,  as  well  as  the  pertinent  values  for 
the  resolvent  kernel  J"1*  itself,  have  been  determined  for  a goodly  number  of 
points  within  the  range  of  x-values  running  from  0 to  4.2,  and  the  results 
have  been  entered  in  the  appended  Table  I. 
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3o.  Determinate  on  of  the  Supersonic  Source  Distribution  Which  Will  Produce 
the  Required  Padial  Components  of  Velocity  Along  the  Duct  Contour 

The  just  determined  values  for  the  resolvent  kernel  / (x)  may  now  be 

used  to  find  the  distribution  of  supersonic  sources  which  are  located  on 
the  duct's  axis  and  the  strength  of  which  is  to  vary  in  the  proper  way  in 
order  that  the  radial  velocity  component  of  the  induced  velocity  along  the 
sought  duct's  surface  connecting  the  end-points  Nj_  and  N2  shall  be  equal  to 
the  values  defined  by  the  relations 


UP-r, 


s/n.  k g . . 

A U<*> 


at  the  various  points  of  the  duct  contour 
running  from  Nj  to  N2,  end  approximated 
by  a seven  term  series,  for  which 
n = 0,1,2, 3,4, 5,6 

where  this  time  the  radial  components 
throughout  the  N]to  N2  interval  are 
approximated  by  a "six  term  series,  for 
which  k = 1,2, 3, 4, 5, 6 


The  way  this  source  distribution  must  behave  is  found  by  first  evaluat- 

o 

ing  the  v(s)  function  defined  as 

Jo 


f^z 


r 


where  it  is  now  true  that 


s'(i)  --  -gii+wi) 7,  --i 


B 


(cos  ue 


sin  tee 


Then  it  follows,  upon  using  the  respective  cosine  and  sine  formulations, 


that  one  must  have 


-*  5/”  7)6  die 
(Q’-e 


cos  K 9 


*-/  /I 


and 


A 

(where  ~~  & has  been  inserted  for  s) 
and  thus  it  may  be  deduced  that 


d a 


Now,  as  is  well  known, 


re1 


**  = 1 faxfreO-Ute')-  cos(r,9')J(ne'\ 


and 


re' 


y=r^"de  - /-?-  [cos (*  e‘)  sm(ice')Jfc.6,\ 


where  the  symbols  $ (t?©')  and  J(n  &')  denote,  respectively,  the  Fresnel 
integrals  of  the  first  and  second  kinds,  with  the  argument  which  is  repre- 
sented by  the  quantity  within  the  parentheses. 

By  use  of  these  Fresnel  integrals,  therefore,  the  expressions  for 
B v(s)  may  be  evaluated;  the  results  are  listed  in  Tables  II  and  III.  Once 
these  values  have  been  obtained,  one  may  easily  find  the  corresponding 
values  for  B w(s),  because  one  has  simply  that 


6 t 
IT^B  ' 


it+/se~ 


0,5198 


8/\r(^y 


This  calculation  has  been  carried  out,  and  the  results  are  given  in 
Tables  IV  and  V. 
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TABLE  II 


EVALUATION  FOR  THE  Bv(s)  FUNCTION  APPEARING  AS  INTERMEDIARY  STAGE 
BETWEEN  CONVERSION  OF  DUCT’S  LOCAL  SLOPE  TO  CORRESPONDING  SOURCE 
INTENSITY  FUNCTION,  FOR  THE  CASE  OF  THE  COSINE  FORMULATION  DESCRIB- 
ING THE  RADIAL  COMPONENTS  OF  VELOCITY;  i.e.,  for  n = C08 

( n l+pB^-e 


X 

0 

1 

2 

3 

4 

5 

6 

0 

OO 

oo 

OO 

oo 

oo 

oo 

0.1 

3.16228 

3.11873 

2.93793 

2.77098 

2.47211 

2.09548 

1.64638 

0.2 

2.23607 

2.11279 

1.74805 

1.16417 

0.39621 

-0.49148 

-1.43143 

0.3 

1.82574 

1.59982 

0.95054 

-0.03474 

-1.16876 

-2.42158 

-3.48744 

0.4 

1.53114 

1.23606 

0.23017 

-1.01217 

-2.40986 

-3.51296 

-4.01278 

0.5 

1.41421 

0.93712 

-0.31084 

-1.87575 

-3.14209 

-3.56185 

-2.83643 

0.6 

1.29099 

0.67213 

-0.82643 

-2.46599 

-3.27643 

-2.58930 

-0.52113 

0.7 

1.19523 

0.42788 

-1.32912 

-2.78369 

-2.77186 

-0.87303 

2.03673 

0.8 

1.11803 

0.19810 

-1.70403 

-2.83747 

-1.76027 

1.11759 

3.86237 

0.9 

1.05409 

-0.02006 

-2.01348 

-2.53900 

-0.42550 

2.34295 

4.28371 

1 

1 

-0.21980 

-2.22178 

-2.00566 

-0.99961 

3.83666 

3.03986 

1.1 

0.9535 

2.2664 

3.8296 

0.7359 

1.2 

0.91287 

-0.53438 

-2.31677 

-0.36349 

3.15361 

2.82064 

-1.86681 

1.3 

0.8771 

3.5265 

1.0781 

-3.7087 

1.4 

0.84515 

-0.93934 

-1.96007 

1.44018 

3.30699 

-0.93155 

-4.12099 

1.5 

0.8165 

2.5228 

-2.6750 

-2.9455 

1.6 

0.79057 

-1.20493 

-1.24469 

2.73111 

1.31045 

-3.68082 

-0.61331 

1.7 

0.7670 

-0.1207 

-3.6854 

1.9367 

1.8 

0.74536 

-1.42374 

-0.30086 

3.02905 

-1  52424 

-2.63821 

3.83854 

1.9 

0.7255 

-2.6455 

-0.9566 

4.2361 

2 

0.70711 

-1.57104 

0.70683 

2.13486 

-3.29222 

1.04308 

3.04435 

2.1 

0.6901 

-3.3648 

2.7742 

0.7171 

2.2 

0.67420 

-1.63106 

1.50254 

0.52036 

-2.85008 

3.77749 

-1.37490 

2.3 

0.6594 

-1.8371 

3.7798 

-3.7428 

2.4 

0.64550 

-1.63821 

2.22994 

-1.32003 

-0.50125 

2.77910 

-4.14910 

2.5 

0.6325 

0.9328 

1.0444 

-2.9513 

2.6 

0.62017 

-1.53518 

2.49345 

-2.62286 

2.22415 

-0.95941 

-0.61799 

2,7 

0.6086 

3.1340 

-2.6947 

1.9598 

2.8 

0.59761 

-1.38599 

2.33838 

-2.91398 

3.50672 

-3.69676 

3.81919 

2.9 

0.5872 

3.2777 

-3.6965 

4.2352 

3 

0.57735 

-1.15797 

1.73393 

-2.08228 

2.43570 

-2.69436 

3.04143 
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TABLE  III 


EVALUATION  FOR  THE  Bv(s)  FUNCTION  APPEARING  AS  INTERMEDIARY  STAGE 
BETWEEN  CONVERSION  OF  DUCT»S  LOCAL  SLOPE  TO  CORRESPONDING  SOURCE 
INTENSITY  FUNCTION,  FOR  THE  CASE  OF  THE  SINE  FORMULATION  DESCRIB- 
ING THE  RADIAL  COMPONENTS  OF  VELOCITY;  i.e.,  FOR  n = 5&..X  S . .. 

TT 


■ 

1 

2 

3 

4 

5 

6 

0 

0 

0 

0 

0 

0 

0 

EES 

0.6283 

1.2474 

1.8451 

2.4103 

2.9325 

3.4032 

B9 

0.8821 

1.7040 

2.4064 

2.9360 

3.2409 

3.2976 

B 

1.0653 

1.9648 

2.5498 

2.6925 

2.3903 

1.6104 

0,4 

1.2049 

2.0760 

2.3317 

1.8817 

0.7416 

-0.8904 

0.5 

1.3117 

2.0497 

1.8515 

0.6633 

-1.2066 

-3.1806 

0,6 

1.3894 

1.9039 

1.1387 

-0.7270 

-2.9036 

-4.3332 

0.7 

4.4408 

1.6728 

0.2758 

-2.0414 

-3.8779 

-3.8940 

0.8 

1.4680 

1.3306 

-0.6296 

-3.0384 

-3.8565 

-2.0226 

0.9 

1,4721 

0.9298 

-1.5044 

-3.5381 

-2.8372 

0.5784 

1 

1.4494 

0.4690 

-2.2490 

-3.4494 

-1.0878 

2.9280 

1.1 

-2.7831 

0.9278 

4.1292 

1.2 

1.3462 

-0.5141 

-3.0640 

-1.6515 

2.6729 

3.7242 

1.3 

-0.2513 

3.6829 

1.8744 

1.4 

1.1829 

-1.4435 

-2.7534 

1.1782 

3.6889 

-0.7086 

1 »5 

2.3907 

2.6915 

-3.0390 

1.6 

0.9409 

-2.1484 

-1.4302 

3.1790 

0.9602 

-4.2282 

1.7 

3.4087 

-1.0412 

-3.8160 

1-8 

0.6574 

-2.5018 

0.4090 

3.0408 

-2.7742 

-1.9572 

1.9 

2.1394 

-3.7743 

0.6366 

2 

0.3317 

-2.4391 

2.0704 

0.8538 

-3,7721 

2.9718 

2.1 

-0.5786 

-2.7676 

4.1646 

2.2 

-0.0143 

-1.9680 

2.9198 

-1.9219 

-1.0303 

3.7578 

2.3 

-2.9389 

0.9760 

1.9044 

2.4 

-0.3635 

-1.1677 

2.6334 

-3.4523 

2.7129 

-0.6852 

2.5 

-3.3735 

3.7168 

-3.0198 

2.6 

-0.7082 

-0.1774 

1.3254 

-2.7180 

3.7185 

■ A .2078 

2.7 

-1.5981 

2.7167 

-3.7968 

2.8 

-1.0207 

0.8331 

-0.5011 

-0.2060 

0.9815 

-1.9440 

2.9 

1.2174 

-1.0220 

0.6474 

3 

-1.2985 

1.6905 

-2.1489 

2.4265 

-2.7567 

2.9862 
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TABLE  IV 


EVALUATION  OF  THE  PART  OF  THE  INTEGRAL  EQUATION,  £0.(36)  OF  PART  I,  WHICH 
MAKES  IT  NONHOMOGEN EOUS,  IN  THE  CASE  OF  THE  COSINE  FORMULATION  FOR  THE 
RADIAL  COMPONENTS  OF  VELOCITY;  i,e., 

TABULATED  VALUES  OF  Bvr(s),  IN  CASE  WHERE  0 n=  — 

e 


* 

0 

1 

2 

3 

4 

5 

6 

0 

OO 

OO 

OO 

OO 

o° 

OO 

OO 

0.1 

1.62351 

1.6011 

1.5340 

1.4226 

1.2692 

1.0758 

0.8453 

0.2 

1.13436 

1.0718 

0.8868 

0.5906 

0.2010 

-0.2493 

-0.7261 

0.3 

0.91524 

0.8020 

0.4765 

-0.0174 

-0.5859 

-3 .2139 

-1.7483 

0.4 

0.79361 

0.6126 

0.1389 

-0.5016 

1.1943 

-1.7410 

-1.9887 

0.5 

0.69310 

0.4593 

-0.1523 

-0.9193 

-1.5399 

-1.7457 

-1.3901 

0.6 

0.62574 

0.3258 

-0.4006 

-1.1953 

-1.5881 

-1.2550 

-0.2526 

0.7 

0.57323 

0.2052 

-0.6374 

-1.3351 

-1.3295 

-0.4187 

0.9768 

0.8 

0.53051 

0.0940 

-0.8085 

-1.3464 

-0.8352 

0.5303 

1.8327 

0.9 

0.49500 

-0.0094 

-0.9455 

-1.1923 

-0.1998 

1.3350 

2.0116 

1 

0.46490 

-0.1022 

-1.0329 

-0.9324 

0.4647 

1.7837 

1.4365 

1.1 

0.43890 

1.0432 

1.7628 

0.3387 

1.2 

0.41618 

-0.2664 

-1.0562 

-0.1680 

1.4377 

1.2859 

-0.8511 

1.3 

0.39610 

1.5926 

0.4869 

-1.6748 

1.4 

0.37804 

-0.4204 

-0.8767 

0.6442 

1.4792 

-0.4167 

-1.8433 

1.5 

0.36195 

1.1184 

-1.1849 

-1.3057 

1.6 

0.34730 

-0.5293 

—0 -5468 

1.1998 

0.5757 

-1.6170 

-0.2697 

1.7 

0.33395 

-0.0526 

-1.6046 

0.8648 

1.8 

0,32170 

-0.6145 

-0.1299 

1.3073 

-0.6579 

-1.1602 

1.6567 

1.9 

0.31051 

-1.1323 

-0.4094 

1.8131 

2 

0.30010 

-0.6667 

0.3000 

0.9273 

-1.3972 

0.4427 

1.2920 

2.1 

0.29046 

-1.4162 

1.1677 

0.3018 

2.2 

0.28148 

-0.6810 

0.6691 

0.2173 

-1.1899 

1.5771 

-0.7828 

2.3 

0.27312 

-0.7609 

1.5656 

-1.5503 

2.4 

0.26530 

-0.6733 

0.9165 

-0.5425 

-0.2060 

1.1422 

-1.7053 

2.5 

0.25793 

0.3804 

0.4259 

-1.2035 

2.6 

0.25098 

-0.6215 

1.0091 

-1.0615 

0.9001 

-0.3883 

-0.2501 

2.7 

0.24441 

1.2586 

-1.0822 

0.7871 

2.8 

0.23821 

-0.5525 

0.9321 

-1.1615 

1.3978 

-1.4735 

1.5223 

2.9 

0.23241 

1.2973 

-1.4631 

1.6763 

3 

0.22684 

-0.4550 

0.7009 

-0*8131 

-0.9766 

-1.0586 

1.1950  ' 
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TABLE  V 


EVALUATION  OF  THE  PART  OF  THE  INTEGRAL  EQUATION, EQ. (36)  OF  PART  I, 
WHICH  MAKES  IT  NONHOMOGEN EOUS, IN  THE  CASE  OF  THE  SINE  FORMULATION 
FOR  THE  RADIAL  COMPONENTS  OF  VELOCITY;  i.e., 

TABULATED  VALUES  OF  B w(s),  IN  CASE  WHERE 


V k 
s 

1 

2 

3 

4 

5 

6 

0 

0 

0 

0 

0 

o 

0 

0.1 

0.3226 

0.6404 

0.9473 

1.2374 

1.5055 

1.7472 

0.2 

0.4475 

0.8644 

1.2207 

1.4893 

1.6440 

1.6728 

0.5 

0.5341 

0.9850 

1.2783 

1.3499 

1.1984 

0.8074 

0.4 

0.5972 

1.0289 

1.1556 

0.9326 

0.3676 

-0.4413 

0.5 

0.6428 

1.0045 

0.9073 

0.3251 

-0.5913 

-1.5587 

0.6 

0.6735 

0.9228 

0.5519 

-0.3524 

-1.4074 

-2.1005 

0.7 

0.6910 

0.8022 

0.1323 

-0.9790 

-1,8597 

-1.8674 

0.8 

0.6966 

0.6314 

-0.2988 

-1.4418 

-1.8300 

-0.9598 

0.9 

0.6913 

0.4367 

-0.7065 

-1.6616 

-1.3324 

0.2716 

1 

0.6738 

0.2180 

-1.0456 

-J .6036 

-0.5057 

1.3612 

1.1 

-1.2811 

0.4271 

1.9008 

1.2 

0.6137 

-0.2344 

-1.3968 

-0.7529 

1.2185 

1.6977 

1.3 

-0.1135 

1.6630 

0.8464 

1.4 

0.5291 

-0.6457 

-1.2317 

0.5270 

1.6502 

-0.3170 

1.5 

1.0597 

1.1931 

-1.3471 

1.6 

0.4134 

-0.9439 

-0.6283 

1.3966 

0.4218 

-1.8576 

1.7 

1.4843 

-0,4534 

-1.6616 

1.8 

0.2838 

-1.0799 

0.1765 

1.3125 

-1.1977 

-0.8448 

1.9 

0.9157 

-1.6154 

0.2725 

2 

0.1408 

-1.0352 

0.8787 

0.3645 

-1.6009 

1.2613 

2.1 

-0.2436 

-1.1650 

1.7530 

2.2 

-0.0060 

-0.8216 

1.2190 

-0.8024 

-0.4302 

1.5689 

2.3 

1.2172 

0.4042 

0.7887 

2.4 

-0.1494 

-0.4799 

1.0822 

-1.4187 

1.1149 

-0.2316 

2.5 

-1.3756 

1.5156 

-1.2514 

2.6 

-0.2866 

-0.0718 

0.5363 

-1.0999 

1.5047 

-1.7027 

2.7 

-0.6419 

1,0911 

-1.5249 

2.8 

-0.4069 

0.3321 

-0.1998 

-0.0321 

0.3913 

-0.7749 

2.9 

0.4818 

-0  .4045 

-0.2562 

3 

-0.5102 

0.6642 

-0.8443 

0.9534 

-1.0831 

1.1733 

m+K 


sinks 
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The  next  step  is  to  insert  the  above-derived  values  for  w(s)  in  Eq.  (37') 

* 

of  Part  I in  order  to  find  the  values  of  B h (s).  The  latter  function  defines 
the  way  in  which  the  pertinent  supersonic  source  distribution  has  to  be  appor- 
tioned along  the  duct's  axis.  The  values  for  this  source  intensity  function 
for  the  specific  case  under  study  have  been  determined  and  the  results  are 
tabulated  in  Tables  VI  and  VII;  in  addition,  they  are  plotted  in  Figs.  1 and 
2,  as  thirteen  individual  curves,  which  correspond  to  the  thirteen  harmonic 
components  used  to  approximate  the  radial  velocity  distribution- 


4.  Computation  of  the  Axial  Components  of  the  Induced  Velocity  Arising  from 
the  Individual  Supersonic  Source  Distributions 

The  axial  components  of  the  velocity  field  induced  at  the  points  on  the 
meridional  contour  running  between  Nj  and  N2  are  now  to  be  calculated  in  the 
case  of  each  one  of  the  thirteen  supersonic  source  distributions  determined 
above.  The  equation  which  defines  these  axial  components  of  the  induced  flow 
in  terms  of  the  source  intensity  distribution  is 


_ 


IL  9 * 


= 6 


h fca/j  


This  formula,  however,  is  not  convenient  to  use  in  performing  the  sctual 
numerical  evaluations  because  of  the  singularity  possessed  by  the  integrand 


at  the  point  for  which  z = s = 0 when  one  has  to  deal  with  the  cosine  formu- 
lation for  the  radial  comoonents,  given  as  = — £os—V Z. , 

(n  7+Jb  *l  & 

Nevertheless  this  difficulty  may  be  easily  side-stepped  by  converting 
the  expression  (4)  into  the  following  rewritten  version,  which  readily  allows 
one  to  carry  out  the  desired  calculation  even  in  this  troublesome  case  of  the 
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TABLE  VI 


ADJUSTED  SUPERSONIC  SOURCE  INTENSITY  FUNCTION,  Bh*(s),  IN  THE  CASE  OF  THE 
COSINE  FORMULATION  FOR  THE  RADIAL  COMPONENTS  OF  VELOCITY,  H n=  — 282*.  e-- 

__  _ _•  1+3B  -y.  © 


\n 

s\ 

0 

1 

2 

3 

4 

5 

6 

0 

©o 

oo 

oo 

oc? 

o© 

oo 

OO 

0.1 

1.52062 

1.49821 

1.43111 

1.31971 

1.16631 

0,97291 

0.74241 

0.2 

0.99547 

0.93428 

0.75335 

0.46372 

0.08287 

-0.35706 

-0.82248 

0.3 

0.75462 

0,64541 

0.33168 

-0.14376 

-0.68915 

-1.29070 

-1.79879 

0.4 

0.60846 

0.44555 

-0.00500 

-0.61103 

-1.26304 

-1.76733 

-1.95627 

0.5 

0.50796 

0.28774 

-0.28611 

-1.00004 

-1,56327 

-1.72854 

-1.33109 

0,6 

0.43366 

0,15418 

-0.51750 

-1.23926 

-1.56460 

-1.18645 

-0.17437 

0.7 

0.37641 

0.03704 

-0.73203 

-1.34101 

-1.26636 

-0.33140 

1.03883 

0.8 

0.33059 

-0.06784 

-0.87716 

-1.31486 

-0.74643 

0.60976 

1.85095 

0.9 

0.29323 

-0,16260 

-0.98621 

-1.12846 

-0 .10303 

1.38339 

1.97589 

1 

0.26223 

-0,24490 

-1.04499 

-0.84450 

0.55141 

1.78694 

1.35690 

1,1 

0.23639 

-0,31453 

-1.03098 

-0,45616 

1.10472 

1.71946 

0.24131 

1,2 

0,21468 

-0,38416 

-1.01698 

-0.06781 

1.46483 

1.20742 

-0.93562 

1.3 

0.19590 

-0.44593 

-0,90520 

0.32681 

1.57912 

0.39041 

-1.72115 

1.4 

0.17971 

-0.50770 

-0,79343 

0.72143 

1.42888 

-0.50765 

-1.34275 

1.5 

0.16563 

-0.54616 

-0,61618 

0.97347 

1.03975 

-1.25061 

-1.26686 

1,6 

0.15334 

-0.58462 

-0.43894 

1.22551 

0.48170 

-1.64514 

-0.21461 

1,7 

0.14244 

-0,61118 

-0.22836 

1.24971 

-0.14727 

-1.59338 

0.90907 

1.8 

0,13287 

-0.63775 

-0,01779 

1.27391 

-0.73895 

-1.11792 

1.66871 

1.9 

0.12432 

-0.64844 

0.18968 

1.06135 

-1.18868 

-0.35212 

1.78444 

2 

0.11672 

-0,65912 

0.39714 

0.84879 

-1.42241 

0.49538 

1.22991 

2.1 

0.10985 

-0,65214 

0.56664 

0.48574 

-1.40946 

1.19981 

0.22543 

2,2 

0.10376 

-0.64517 

-0.73615 

0.12268 

-1.15585 

1.57730 

-0.84934 

2.3 

0.09799 

-0,62890 

0,84034 

-0.24918 

-0.70822 

1.53192 

-1.58862 

2.4 

0.09313 

-0,61264 

0.94452 

-0.62104 

-0.14606 

1.08197 

-1.70691 

2.5 

0.08853 

-0,57633 

0.97027 

-0.36086 

0.45517 

0.35266 

-1.17471 

2.6 

0.08436 

-0.54002 

0.99602 

-1.10067 

0.93889 

-0,45797 

-0.20785 

2.7 

0.08047 

-0.49776 

0,93918 

-1.12754 

1,27367 

-1.13335 

0,82176 

2.8 

0.07700 

-0.45549 

0.88234 

-1.15441 

1.38586 

-1.49650 

1.53222 

2,9 

0.07381 

-0,40133 

G. 75340 

-0.96468 

1.25973 

-1.45611 

1.65422 

3 

0.07086 

-0.34717 

0.62447 

-0.77496 

0.91915 

-1.02753 

— 

1.14612 
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TABLE  VII 


ADJUSTED  SUPERSONIC  SOURCE  INTENSITY  FUNCTION, Bh*(s),  IN  THE  CASE 
OF  THE  SINE  FORMULATION  FOR  THE  RADIAL  COMPONENTS  OF  VELOCITY, 

fl  = sin  kQ 
(tm-k 


\ k 
s \ 

1 

2 

3 

" " ’ 

4 

5 

6 

0 

0 

0 

0 

0 

0 

0 

0.1 

0.31719 

0.69267 

0.93145 

1.21665 

1.48025 

\ 

1.71791 

0.2 

0.43018 

0.83041 

1.17137 

1.42659 

1.57045 

1.59129 

0.3 

0.50274 

0.92467 

1.19368 

1.24766 

1.08620 

0.69364 

0.4 

0.55083 

0.94221 

1.04002 

0.80308 

0.23990 

-0.55202 

0.5 

0.58117 

0.89356 

0.76911 

0.18606 

-0.70648 

-1.33160 

0.6 

0.59685 

0.79129 

0.40189 

-0.43145 

-1.48548 

-2.11491 

0.7 

0.59991 

0.65475 

-0.01751 

-1.03047 

-1.88605 

-1.82477 

0.8 

0.59196 

0.47341 

-0.43661 

-1.50295 

-1.80338 

-0.88105 

0.9 

0.57427 

0.27399 

-0.82193 

-1.67663 

-1.26470 

0.35403 

1 

0.54555 

0.05645 

-1.13056 

-1.57426 

-0.41817 

1.41593 

1.1 

0.50754 

-0.15748 

-1.26650 

-1.21610 

0.50968 

1.90886 

1.2 

0.46952 

-0.37140 

-1.40244 

-0.66603 

1.27443 

1.65886 

1.3 

0.42130 

-0.55862 

-1.28785 

-0.02111 

1.67904 

0.77812 

1.4 

0.37308 

-0.74584 

-1.17325 

0.60870 

1.62433 

-0.38715 

1.5 

0.31249 

-0.87086 

-0.85320 

1.11743 

1.13457 

-1.39203 

1.6 

0.25190 

-0.99588 

-0.53314 

1.42192 

0.34793 

-1.86091 

1.7 

0.18757 

-1.03814 

-0.13080 

1.47468 

-0.52206 

-1.62298 

1.8 

0.12324 

-1.08041 

0.27154 

1.27151 

-1.24292 

-0.76512 

1.3 

0.05518 

-1.03472 

0.6C622 

0.85203 

-1.62583 

0.34052 

2 

-0.01288 

-0.98903 

0.94089 

0.29029 

-1.57469 

1.30739 

2.1 

-0.08018 

-0.86465 

1.08583 

-0.31489 

-1.10986 

1.76350 

2.2 

-0.14747 

-0.74027 

1.23076 

-0.35828 

-0.36132 

1.54249 

2.3 

-0.21084 

-0.55980 

1.13836 

-1.24875 

0.46962 

0.73831 

2.4 

-0.27420 

-0.37952 

1.04595 

-1.42108 

1.16055 

-0.33466 

2.5 

-0.33264 

-0.17433 

0.76438 

-1.34919 

1.53135 

-1.26585 

2.6 

-0.39107 

0.03066 

0.48281 

-1.04986 

1.48863 

-1.70523 

2.7 

-0.43972 

0.22556 

0.10910 

-0.57711 

1.04967 

-1.49459 

2.8 

-0.48037 

0.42046 

-0.26461 

-0.01362 

0.33731 

-0.72308 

2.9 

-0.52770 

0.57330 

-0.57409 

0.54344 

-0.45553 

0.31071 

3 

-0.56703 

0.72615 

-0.88356 

0.99751 

-1.11707 

1.21142 
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HARMONIC  COMPONENTS  OF  THE  ADJUSTED  SUPERSONIC  SOURCE 

\ 

i INTENSITY  FUNCTION,  Bh'ts),  IN  THE  CASE  OF  THE  COSINE  FORMULATION 

cos  n 0 

FOR  THE  RADIAL  COMPONENTS  OF  VELOCITY \ 1?,  * j 

1+JDB  -ijfr  0 


HARMONIC  COMPONENTS  OF  THE  ADJUSTED  SUPERSONIC  SOURCE 
INTENSITY  FUNCTION,  Bh*(s),  IN  THE  CASE  OF  THE  SINE  FORMULATION 
FOR  THE  RADIAL  COMPONENTS  OF  VELOCITY,  V-J  ~ t 


I REDUCED  DISTANCE  ALONG  DUCT  AXIS 

FIGURE  2 


cosine  formulation: 


B 


Mjl-.sPJlj 

s * l 


Uo  9k 


, YT 


(e) 


,** 


where  the  h (s)  function  is  defined  as 


Bf)  (■*> ' ~p7%2 


f • 

YTYee  ' 


/-/SB  + /-/3B 


(6) 


"'Rf 

This  newly  introduced  h'  (s)  function  is  now  continuous  at  the  point  for 
which  z = s = 0,  and  thus  the  indicated  quadrature  involving  this  function 
that  appears  on  the  right  hand  side  of  Eq.  (5)  may  be  carried  out  numerically 
without  further  difficulty*.  The  second  integral  appearing  on  the  right  hand 
side  of  FOo  (5)  may  be  evaluated,  of  course,  in  closed  form,  by  use  of  the 
complete  elliptic  integral  of  the  first  kind*  Thus  the  evaluation  of  the 
axisi  velocity  components  in  the  case  where  the  cosine  formulation  is  being 
treated  may  be  carried  out  by  use  of  the  relations: 


B_ 

Of# 


(7) 


where 


f„{j)  = P ej/\±)  + ==L 

J /Z666 


03  96 


6667^  t-2.66^7  i [/ 4*667} + 


IVs) 


in  which  B H (s)  has  the  specific  fora 


ni*t,  x &h  *U>)  - 0.5 i 96  • yip 

oh  - j ■ ■ ■ i i'"" 

J/.66(o7  y - +■  2.G>ic&  7 


These  evaluations  for  Fn(z)  have  been  performed  in  the  case  of  interest 
and  the  results  are  presented  in  Table  VIII,  and,  in  addition,  the  results 
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TABLE  VIII 


+B  5 n 

AXIAL  COMPONENTS  OF  INDUCED  VELOCITIES  AT  THE  DUCT'S  SURFACE, Fn=  U„  8x  * 
CALCULATED  BY  MEANS  OF  Eq=(47)  OF  PART  I, AND  COMPARISON  OF  THESE  VALUES  WITH 
THOSE  OBTAINED  FROM  THF  PROCEDURE  EMPLOYING  ACKERET'S  TWO-DIMENSIONAL  FORMULA , 
IN  WHICH  CASE  Fn*=  hfjsB^d  THE  RADIAL  COMPONENTS  ARE  DESCRIBED  BY  MEANS 

OF  THE  COSINE  FORMULATlBf  FOR  WHICH  0*  = C°fAt 


Fo*=  1 
Fo  = 1 

0.94356 

0=9756 

0.88898 

0 = 9524 

0.83^82 

0=9302 

0.79525 

0.9091 

0.75467 

0.8889 

0=71751 

0=8696 

0.68345 

0=8511 

0=65207 

0.8333 

0.62311 

0=8163 

0=59631 

0=8000 


0.54871 

0.7692 


0=50774 

0.7407 


0.47222 

0.7143 


0.44121 

0.6896 


FV1 

Fi  =1 

0.93713 

0.9703 

0 = 86747 

0.9316 

0.79513 

0.8847 

0.72031 

0.8305 

0,64338 

0.7698 

0.56478 

0.7035 

0.48512 

0.6325 

0.40494 

0=5576 

0=32489 

0.4798 

0.24646 

0.4000 


.09080 

.2377 


-0.05417 

0.0774 


-0.18277 

-0.0747 


-0.29399 

-0.2131 


3N  FOR  WHICH  nH=-r 
_I_ 12 u, 

2 3 

f^ti  f5  =i 

E2  =1  Fj*=1 

0.91788  0.88592 
0.9543  0.9279 

0=80372  0.70094 
0.8700  0=7705 

0.66546  0.46485 
0.7526  0.5468 

0.50889  0.20682 
0.6083  0.2809 

0.34215  -0.06030 
0.4444  0 

0.17352  -0.30426 
0.2687  -0.2687 

0.00193  -0.50438 
0=0089  -0.5002 
-0=15764  -0.64839 
-0.0871  -0.6742 

-0.30342  -0.71705 
-0.2523  -0.7764 

-0.42971  -0.71236 
-0.4000  -0.8000 


-0.61128  -0.49212 
-0.6223  -0.6223 


-1.66604  -0.10096 
-0.7246  -0.2289 


7~ 

3 

4 

5 

6 

FS  =1 
f3*=i 

F4  =1 
F4*=l 

•V1 

F5  =1 

F6  =1 

0.88592 

0.84191 

0.78642 

0.72029 

0.9279 

0.8913 

0.8449 

0.7893 

0.70094 

0.56449 

0.40387 

0=22914 

0=7705 

0.6372 

0.4762 

0.2943 

0.46485 

0.22256 

-0.05144 

-0.30866 

0.5468 

0.2875 

0 

-0=2875 

0.20682 

-0.13645 

-0.46055 

-0.68734 

0.2809 

-0.0950 

-0.4545 

-0=7355 

-0.60061 

-0.6987 


-0.43570 

-0.5579 


0.51108 

0.2207 


0.59482 

0.5579 


-0.45125 

-0.4444 

-0.67105 

-0.7035 

-0.75908 

-0.8325 

-0.70794 

-0.8151 

-0=53357 

-0.6604 

-0.27088 

-0.4000 

0.03234 

-0=0820 

0.32351 

0.2377 

0.55552 

0.5050 

0.69151 

0.6767 

0.71102 

0.7273 

0.61422 

0.6525 

0=42112 

0.4696 

0.16725 

0.2131 

-0=10209 

-0=0709 


-0.72034 

-0.7698 

-0.77030 

-0.8696 

-0.60690 

-0.7370 

-0.28196 

-0.4167 

-0.11292 

0 

0=47168 

0.4000 

0=70195 

0.6792 

0.74844 

0=7692' 

0.60510 

0.6536 

0.31656 

0.3704 

-0.03694 

0 

-0.36090 

-0.3571 

-0.57185 

-0.6077 

-0.61889 

-0.6897 

-0=49630 

-0.5871 


-0.77822 

-0.8888 

-0.55803 

-0.7035 

-0.12557 

-0.2630 

-0=34837 

0.2575 

0.68888 

0.6604 

0.77083 

0.8000 

0.57176 

0.6345 

0.17692 

0.2377 

-0.25679 

-0.2332 

-0.56877 

-0.5993 

-0.64864 

-0.7273 

-0.47375 

-0=5779 

-0.11764 

-0.2168 

0.27827 

0.2131 

0.56262 

0.5485 


TABLE  VIII  (Cont.) 


\ 

0 

1 

2 

3 

4 

5 

6 

2 

0.41382 

0.6667 

-0.38564 

-0.3333 

-0.20577 

-0.3333 

0.65169 

0.6667 

-0.34021 

-0.3335 

-0.24208 

-0.3333 

0.63398 

0.6667 

2.1 

-0.50726 

-0.5303 

0.07488 

0 

0.47585 

0.5305 

2.2 

0.38990 

0.6452 

-0.45277 

-0.4317 

0.04882 

-0.0674 

0.47582 

0.5219 

-0.58040 

-0.6311 

0.36602 

0.3226 

0.14930 

0.1994 

2.3 

-0.54948 

-0.6210 

0.55682 

0.5497 

-0.21735 

-0.1962 

2.4 

0.36847 

0.6250 

-0.49865 

-0.5056 

0.28317 

0.1931 

0.14351 

0.1931 

-0.42308 

-0.5056 

0.60008 

0.6250 

-0.48567 

-0.5056 

2.5 

-0.22528 

-0.3077 

0.48863 

0.5329 

-0.55314 

-0.6154 

2.6 

0.34907 

0.6061 

-0.51800 

-0.5537 

0.46201 

0.4055 

-0.21180 

-0.1872 

0.00853 

-0.0634 

0.25721 

0.3030 

-0.40556 

-0.4903 

2.7 

0.23860 

0.1845 

-0.03202 

0 

-0.1Q716 

-0.1845 

2.8 

0.33210 

0.5882 

-0.51485 

-0.5754 

0.55910 

0.5374 

-0.46033 

-0.4759 

0.42392 

0.3936 

-0.29910 

-0.2941 

0.23067 

0.1818 

2.9 

0.53490 

0.5296 

-0.47539 

-0.5020 

0.47823 

0.4690 

3 

0.31661 

0.5714 

-0.48915 

-0.5714 

0.56210 

0.5714 

-0.51861 

-0.5714 

0.55443 

0.5714 

-0.51702 

-0.5714 

0.54413 

0.5714 

» 


AXIAL  COMPONENTS  OF  INDUCED  VELOCITY  AT  THE  DUCT'S  SURFACE, 
F„*-&  .CALCULATED  BY  MEANS  OF  EQ.  (47)  OF  PART  I AND 

COMPARISON  OF  THIS  CURVE  WITH  THE  ONE  OBTAINED  FROM  THE 
PROCEDURE  EMPLOYING  ACKEREfS  TWO-DIMENSIONAL  FORMULA,  IN 
WHICH  CASE  Fh*  HMEV  THE  RADIAL  COMPONENTS 

ARE  DESCRIBED  BY  MEANS  OF  THE  COSINE  FORMULATION , FOR 
WH,CH 


FIGURE  3 b 
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AXIAL  COMPONENTS  OF  INDUCED  VELOCITY  AT  THE  DUCT'S  SURFACE, 
F**~U  ‘ CALCULAT£D  BY  MEANS  OF  EO.  (47)  OF  PART  I AND 

COMPARISON  OF  THIS  CURVE  WITH  THE  ONE  OBTAINED  FROM  THE 
PROCEDURE  EMPLOYING  ACKERET'S  TWO-DIMENSIONAL  FORMULA,  IN 
WHICH  CASE  Fh-  WHEN  THE  RADIAL  COMPONENTS 

ARE  DESCRIBED  BY  MEANS  OF  THE  COSINE  FORMULATION , FOR 

u/lhs*u  n — COS 
wh/ch 


FIGURE  3d 


AXIAL  COMPONENTS  OF  INDUCED  VELOCITY  AT  THE  DUCfS  SURFACE, 
Fn‘-&  -^LL , CALCULATED  BY  MEANS  OF  EQ.  (47)  OF  PART  I AND 
COMPARISON  OF  THIS  CURVE  WITH  THE  ONE  OBTAINED  FROM  THE 
PROCEDURE  EMPLOYING  ACKERET'S  TWO-DIMENSIONAL  FORMULA,  IN 
WHICH  CASE  Fh-  WHEN  THE  MDtAL  COMPONENTS 

ARE  DESCRIBED  BY  MEANS  OF  THE  COSINE  FORMULATION,  FOR 

mCH  ■ V 


FIGURE  3e 


AXIAL  COMPONENTS  OF  INDUCED  VELOCITY  AT  THE  DUCT'S  SURFACE, 
F^-rr  4^-  .CALCULATED  BY  MEANS  OF  EO.  (47)  OF  PART  I AND 

l/ag  OX 

COMPARISON  OF  THIS  CURVE  WITH  THE  ONE  OBTAINED  FROM  THE 
PROCEDURE  EMPLOYING  ACKEREfS  TWO-DIMENSIONAL  FORMULA,  IN 
WHICH  CASE  Fn*  m!EN  THE  RADIAL  COMPONENTS 

ARE  DESCRIBED  BY  MEANS  OF  THE  COSINE  FORMULATION , FOR 

WH,CH  V / ~7b%7^ 


0 0.2  0.4  0.6  06  1.0  1.2  1.4  1.6  1.8  2 JO  2 2 2.4  2.6  2.8 

Z’tL-I  reduced  DISTANCE  along  duct  axis 


FIGURE  3g 
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are  displayed  graphically  in  Fig.  3,  for  the  seven  harmonic  components 


being  employed  throughout  these  calculations. 


The  values  of  r 


IL 


d X 


have  also  been  computed  on  the  basis  of 


Ackeret's  two-dimensional  formula  for  each  one  of  the  radial  components,  Kjj,  • 

These  values  ere  simply  given  as  f-  * - — — , and  these  two- 

» l-h/88 

dimensional  results  have  been  entered  into  Table  VIII  likewise,  and,  in  addi- 
tion, they  have  also  been  plotted  for  comparative  purposes  in  Fig.  3. 

In  the  corresponding  calculations  which  have  to  be  made  when  it  comes  to 
use  of  the  sine  formulations  for  the  radial  comocnents,  no  difficulty  such 

as  encountered  above  will  arise,  because  in  the  process  of  carrying  out  the 

s/n  K 0 

numerical  evaluations  in  the  case  where  Y)  = - — ~ - it  will  be  found 


that  the  integral  that  defines 


, in  the  expression  for 


has  no  singularity,  in  this  case,  at  the  point  where  z = s = 0„  The 

straight-forward  numerical  computations  have  been  easily  accomplished  for 

the  Fjc(z)  functions,  consecuently,  and  the  results  are  presented  in  Table  IX, 

while  the  corresponding  graphical  plots  are  given  in  Fig.  4.  The  values  of 
■ft 

Fk  , which  are  defined  according  to  the  relation 


c*  - sm  < & 

' !+/S8  e 

corresponding  to  the  two-dimensional  interpretation  of  the  problem  that 
makes  use  of  Ackeret's  formula  for  the  pressures,  have  also  been  computed, 
and  these  results  have  likewise  been  entered  into  Table  IX  and  plotted  up 
on  the  graphs  of  Fig.  4,  for  comparison  with  the  more  exact  three-dimensional 
values. 
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TABLE  IX 


AXIAL  COMPONENTS  OF  INDUCED  VELOCITIES  AT  THE  DUCT'S  SURFACE, Ffc*  +&-  , 

CALCULATED  BY  MEANS  OF  EQ.(47)  OF  PART  I, AND  COMPARISON  OF  THESE  VALUES  WITH 


THOSE  OBTAINED  FRQMTHE  PROCEDURE  EMPLOYING  ACKERET'S  TWO-DIMENSIONAL  FORMULA 
IN  WHICH  CASE  Fk*=  , WHEN  THE  RADIAL  COMPONENTS  ARE  DESCRIBED  BY 


l*f®4  0 1 

MEANS  OF  THE  SINE  FORMULATION  FOR  WHICH 


_ ain  k 9 


X 

1 

2 

3 

4 

5 

6 

0 

Fi  = 0 
Fi*=  0 

F2  = 0 
F2*=  0 

F3  = 0 
F3*=  0 

F4  = 0 
F4*=  0 

F5  = 0 
F5*=  0 

F6  = 0 
F6*=  0 

0.1 

0.09101 

0.1020 

0.13066 

0.2028 

0.26724 

0.3015 

0.34907 

0.3968 

0.42470 

0.4878 

0.49289 

0.5754 

0.2 

0.13084 

0.1980 

0.35233 

0.3874 

0.50502 

0.5598 

0.63032 

0.7078 

0.72005 

0.8248 

0.77009 

0.9058 

0.3 

0.26033 

0.2875 

0.49191 

0.5468 

0.66840 

0.7526 

0.76681 

0.3847 

0.78164 

0.9302 

0.70584 

0.8847 

0.4 

0.33038 

0.3698 

0.59767 

0.6756 

0.74732 

0.8646 

0.75269 

0.9041 

0.61234 

0.7873 

0.34851 

0.5343 

0.5 

0.39137 

0.4444 

0.66697 

0.7698 

0.74423 

0.8889 

0.59885 

0.7698 

0.27015 

0.4444 

-0.15203 

0 

0.6 

0.44340 

0.5111 

0.69938 

0.8270 

0.66220 

0.8270 

0.33841 

0.5111 

-0.14594 

0 

-0.59830 

-0.5111 

0.48657 

0.5695 

0.69835 

0.8464 

0.51221 

0.6885 

0.01980 

0.1769 

-0.52274 

-0.4255 

-0.88395 

-0.8094 

0.8 

0.52105 

0.6193 

0.66228 

0.8288 

0.31335 

0.4898 

-0.30007 

-0.1732 

-0.76216 

-0.7217 

-0.75262 

-0."925 

0.9 

0.54700 

0.6604 

0.53652 

0.7764 

0,08539 

0.2523 

-0.56649 

-0.4798 

-0.80606 

-0.8163 

-0.42121 

-0.4798 

1 

0.56412 

0.6928 

0.50405 

0.6228 

-0.14830 

0 

-0.73607 

-0.6928 

-0.64979 

-0.6928 

0.03864 

0 

1.1 

-0.78330 

-0.7800 

-0.34093 

-0.3921 

0.45175 

0.4610 

1.2 

0.57443 

0.7316 

0.25538 

0.4521 

-0.54885 

-0.4521 

-0.70417 

-0.7316 

0.03494 

0 

0.66659 

0.7316 

1.3 

-0.51640 
—0  *5609 

0.37799 

0.3773 

0.61183 

0.7178 

1.4 

0.55860 

0.7367 

-0.02457 

0.1540 

-0.73017 

-0.7045 

-0.25508 

-0.3013 

0.60048 

0.6415 

0.31922 

0.4354 

1.5 

0.03359 

0 

0.64935 

0.7273 

-0.09252 

0 

1.6 

0.51381 

0.7104 

-0.29319 

-0.1485 

-0.65010 

-0.6793 

0.30060 

0.2905 

0.51862 

0.6186 

-0.46496 

-0.4198 

1.7 

0.50223 

0.5215 

0.24936 

0.3509 

-0.65924 

-0.6674 

» 


I 
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TABLE  IX  (Cont.) 


X 

1 

2 

3 

4 

5 

6 

1-8 

0-44713 

0.6559 

-0.50572 

-0.4054 

-0.35261 

-0.4054 

0.60715 

0-6559 

-0.08265 

0 

-0.60458 

-0.6559 

1.9 

0.60259 

0.6742 

-0.38622 

-0.3390 

-0.33436 

-0.3985 

2 

0.36303 

0.5775 

-0.62876 

-0.5773 

0.04410 

0 

0.49100 

0.5773 

-0.58494 

-0.5775 

-0-04483 

0 

2-1 

0.29433 

0.3854 

-0.62984 

-0.6557 

0.38827 

0.3854 

2.2 

0.26438 

0-4794 

-0.65014 

-0.6416 

0.38613 

0.3792 

0.05173 

0.1341 

-0.51123 

-0.5587 

0.57098 

0.6136 

2.3 

-0.19558 

-0.1320 

-0.26648 

-0.3175 

0.52802 

0.6038 

2.4 

0-15833 

0.3674 

-0.56822 

-0.5944 

0.55742 

0.5944 

-0.40527 

-0.3674 

0.03639 

0 

0.28297 

0.3674 

2.5 

-0.54285 

-0.5329 

0.31626 

0.3077 

-0.06570 

0 

2.6 

0.04837 

0-2465 

-0.40303 

-0.4504 

0.50863 

0.5764 

-0.58884 

-0.6027 

0.49926 

0.5249 

-0.38338 

-0.3562 

2.7 

-0.53952 

-0-5678 

0.54374 

0.5970 

-0.55262 

-0.5678 

2.8 

-0.05999 

0.1223 

-0.18585 

-0.2392 

0.26589 

0.3457 

-0.39336 

-0.4371 

0.45892 

0.5094 

-0.51525 

-0.5594 

2.9 

-0.19295 

-0.2358 

0.21316 

0.2898 

-0.28595 

-0.3407 

3 

-0.16295 

0 

0.04540 

0 

-0.06629 

0 

0.03206 

0 

-0.05697 

0 

0.03871 

0 

AXIAL  COMPONENTS  OF  INDUCED  VELOCITY  AT  THE  DUCT'S  SURFACE, 
F*‘TL  , CALCULATED  BY  MEANS  OF  EO.  (47)  OF  PART  I AND 

COMPARISON  OF  THIS  CURVE  WITH  THE  ONE  OBTAINED  FROM  THE 
PROCEDURE  EMPLOYING  ACKERETS  TWO-DIMENSIONAL  FORMULA,  IN 
WHICH  CASE  Fk*  WHEN  ™E  RADIAL  COMPONENTS 

ARE  DESCRIBED  BY  MEANS  OF  THE  SINE  FORMULATION,  FOR 

WHICH  VwtbjNfysi 
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FIGURE  4 b 
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AXIAL  COMPONENTS  OF  INDUCED  VELOCITY  AT  THE  DUCTS  SURFACE, 
F*'Z l , CALCULATED  BY  MEANS  OF  EQ.  (47)  OF  PART  I AND 

COMPARISON  OF  THIS  CURVE  WITH  THE  ONE  OBTAINED  FROM  THE 
PROCEDURE  EMPLOYING  ACKERETS  TWO-DIMENSIONAL  FORMULA,  IN 
WHICH  CASE  Fk-  /£)$•  WHEN  THE  RADIAL  COMPONENTS 

ARE  DESCRIBED  BY  MEANS  OF  THE  SINE  FORMULATION,  FOR 

WH,CH  Kirfftfos* 
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AXIAL  COMPONENTS  OF  INDUCED  VELOCITY  AT  THE  DUCT'S  SURFACE, 
Fk.J.  , CALCULATED  BY  MEANS  OF  EO.  (47)  OF  PART  I AND 
COMPARISON  OF  THIS  CURVE  WITH  THE  ONE  OBTAINED  FROM  THE 
PROCEDURE  EMPLOYING  ACKERET'S  TWO-DIMENSIONAL  FORMULA,  IN 
WHICH  CASE  Fk‘  WHEN  ™E  RADIAL  COMPONENTS 

ARE  DESCRIBED  BY  MEANS  OF  THE  SINE  FORMULATION,  FOR 
WH,CH 


It  may  immediately  be  seen  upon  comparison  of  the  Fn  values  with  the 
Fn  values,  and  of  the  Fjj  values  with  the  F^  values  that,  throughout  the 
whole  z-range  extending  from  0 to  5.0,  these  values  do  not  differ  from  one 
another  to  any  appreciable  degree  except  for  the  smallest  values  of  n and  k. 

It  is  worthwhile  to  point  out  that,  in  the  case  where  0 = 0,  the  corres- 
ponding  F and  F values  will  become  equal  as  n and  k become  infinite.  On  the 
basis  of  this  observation  one  may  cautiously  presume  that  the  same  result 
will  hold  true  even  when  the  value  of  p is  not  zero,  though  still  of  small  size 
but  in  the  cases  which  come  up  in  most  practical  applications  it  will  be  ade- 
quate, in  general  to  consider  merely  such  suitably  small  values  of  p. 

To  make  the  above  statement  clear  in  the  case  where  p = 0,  one  may  pro- 
ceed to  determine  formally  the  field  of  flow  which  is  external  to  an  infinitely 
long  circular  cylinder,  by  use  of  the  three-dimensional  approach  already 
illustrated.  In  this  case  the  axis  of  the  cylinder  is  taken  to  be  along  the 
x-axis  and  its  radius  is  to  be  equal  to  B Rq,  so  that  the  radial  coordinate 

is  defined  as  r = ^ ° 

d no 

The  n'th  one  of  the  harmonic  components  defining  the  radial  velocity 
produced  by  action  of  the  corresponding  supersonic  source  distribution,  which 
is  arrayed  along  the  axis  of  the  cylinder  in  question,  has  the  magnitude 
denoted  by  cos  (nz)  at  the  surface  (i.e.,  at  all  points  ranged  along 

the  generatrix  of  the  cylinder  lying  in  any  one  of  the  arbitrary  general 
meridional  planes  passing  through  its  axis).  The  velocity  potential  which 
describes  such  a flow  produced  by  the  proper  distribution  of  supersonic 
sources  placed  along  the  axis  is  given  by  means  of  the  formula 


where  JQ  and  Y0  are  the  cylindrical  Bessel  functions  of  the  first  and 
second  kinds,  and  depending  on  the  parameter  nBr,  while  the  symbols  A and 
C stand  for  constants  which  are  determined  by  the  boundary  condition.  This 
boundary  condition  to  be  satisfied  is.  merely: 


It  is  found,  upon  applying  this  condition,  that 


tlA, ) 
' y;(s„) 


A 


and 


tine) 

Y'  (71 6)  + j0'  (/)  6) 


It  follows,  consequently,  that  the  proper  expression  for  the  axial 
component  of  the  auxiliary  potential  is  given  by 


B 


Y.Y.'*JjJ 

Yt+Jt 


s/a 


,J_  Yo  Jg  Y, 


B 


t 


COS 


at  the  surface  of  the  cylinder* 

But  if  the  limiting  values  that  these  coefficients  of  sin(nz)  and 

cos  (nz)  take  on  as  n — > are  examined  it  is  found  that 

cos  % cos  y,  + s/%  % s/n 

-k~  rr  71 8 ~ 

/ 

ir  ?i  Q 

- - cos(Yj-f')  = 0 


Mm  X-YL^tXA a = 

Ye'  + J'  9 *—?>*> 
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because,  as  7)  — ? cfi 


it  is  true  that 


y°M~=7&3 


^8)-w$r 


and 


W-7&& 


herein  % = “ ***  ^ = * & ~(/^)  “ff 

In  like  manner  it  may  also  be  seen  that 

= SIJI  _ %)  = / 

T0  + J0 

and  thus,  for  sufficiently  large  values  of  n,  it  follows  that 

3j^_)  ~_L 

3*  Li 


COS  7! 


b 7 

which  is  identical  with  the  result  that  one  comes  to  in  the  two-dimensional 
treatment  through  use  of  the  Ackeret  formula  for  relating  the  pressure  to 
the  geometrical  characteristics  of  the  duct. 
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50  Evaluation  of  the  Constants  Am  and  Am+  k 


* * 

Having  obtained  the  values  of  Fjj,  Fjj  , Fk,  and  Fj^  , the  next  step  in  the 
procedure  described  in  Part  I is  to  determine  the  quantities  £„  and  Em+k* 

In  the  particular  example  under  consideration  the  appropriate  formulae  for 

carrying  out  the  determinations  are  ^ 

£,  s 2 « 2.19325 

“r1 


and 


W 1 £ 2.19325  [%-£)  y)<i  y 
it*  0 


(10) 


It  should  be  noted  that  in  the  numerical  evaluation  of  the  e-values  one 
needs  to  coraoute  the  integral  f C0S  as  well  as  the  integral 

-J 

f — 13. Q ^.0..  an(i  these  comoutations  can,  if  desired,  be  carried  out  by 

means  of  the  procedure  described  in  Part  I«  As  a practical  way  of  handling 
this  computational  chore,  however,  it  will  be  more  convenient  to  proceed  by 
evaluation  of  the  integrals  written  on  the  right  hand  side  of  the  Eqs.(lO). 

As  a matter  of  fact,  the  e-values  actually  arise  as  the  difference  between 
two  quantities  which  differ  but  slightly  from  one  another.  If  the  integrals 
of  the  Fn  functions  were  to  be  carried  out  by  the  approximate  method  (and  this 
aoproximation  will  naturally  become  less  accurate  as  n grows  to  large  values) 
and  if  the  integrals  of  the  Fn*  functions  were  to  be  obtained  exactly,  then 
the  resulting  error  which  would  be  introduced  into  the  corresponding  values 
for  e would  actually  be  larger  than  the  errors  which  are  made  in  carrying  out 
the  integrations  on  the  Fn-Fn  differences  by  just  numerical  means. 
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Such  computations  have  been  made  and  the  results  are  as  follows  for  the 
specific  example  under  consideration: 


E0  = 1.34424 

£7 

= 0.67938 

£X  = 0.25518 

e8 

= 0.27054 

£2  = 0.17318 

es 

= 0.23286 

£3  £ 0.16114 

£10 

= 0.19772 

£4  = 0.15921 

eU 

* 0.18 

£5  = 0.15728 

£12 

= 0.17 

= 0.15535 

Now  the  coefficients  y//  . which  are  employed  in  Eqs.  (59')  and  (59") 
of  Part  I have  the  values  which  are  displayed  in  the  appended  Table  X,  so 
that  the  equations  which  are  used  in  determining  the  sought  An  values  take 
the  form: 

E0  = 1.34424 


9,0 


A -^3.*  A 
* - 2 


9.9 


^9,9 


^•,6 

^9.9 


° ^"tll  2 It 


A -diiUL 

* -*** 


= 3. 6577Ao-0.82166A2-0.15299A4-0.06183A6 +(1.62114)  (0.67938) 
■$.06004)  (0.23286)  + (0. 06485)  (0.18) . 
ea  = 0.17318 


- 37  - 


ARRAY  OF  COEFFICIENTS,>fa,i,  APPEARING  IN  EQUATIONS  DETERMINING  e*s,  ACCORDING  TO  EC.(S9i 

OF  PART  I,  AND  THEMSELVES  DEFINED  IN  EQ.(58)  OF  PART  I 


0.11755  0 0.19753  0 1.19008  000000  2.46740 


*>*  (a,  *'  su„  ^ 'V 

t/*z  „ ('2f1'"  ~7TlL 'mJL~ ~ ^2t'" 

3-i,t\^lUlll  y^nttl  ^iltH  y'7/>//  ^tltfl 


~ -0 . 82164A0+1 . 52611A2-0 . 54543A4-0 . 10536Ag+ (0 . 18012) (0.67933) 
+(0.58361) (0.23286)+(0. 09190) (0.13) . 
e4=  0.15921 

= -0.15301  Ao-0.54543A2+1.63079A4-0.42823A6+(0. 00721) (0.67938) 
+(0.29776) (0.23286)+(0. 50035) (0.18) . 

Eg=  0.15535 

= -0.06183Ao-0.10535A2-0.42824A4+2.19959Ag+(0.00132) (0.67938) 
+(0.02001) (0.23286)+(0. 33495) (0.18) . 
tX-  0.25518 

= 1.14810Ai-0.52111A3-0.17392A5+(0. 72051) (0.27054) 

(0.11528) (0.19772)+(0. 04764) (0.17) . 

E5=  0.16114 

= -0 . 6211CAx+l . 59420A5-0 . 52294Ag+(0 .25938) (0.27054) 

+(0.52935) (0.19772)+(0. 08006) (0.17) . 
e5=  0.15728 

= -0.17391Ai-0.52^93A3+1.65986A5+(0. 01470) (0.27054) 

+(0.32022) (0.19772)+(0. 048232) (0.17) . 

These  ecuations  then  reduce  down  into  the  simple  forms 

3.36577Ao-0.82166A2-0.15299A4-0.06183Ag  = 0.21722 

0 . 82164A0-1 . 52611A2+0 . 54543A4+0 . 10536Ag  = 0.10163 

0.15301Ao+0.54543A2-1.63079A4+0.42823A6  = 0.00509 

0 .06183Ao+0 .10535A2+0 . 42824A4-2 .13959 Aq  = -0.08950 
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1 »14810Ax-O  . 62111A3-0 . 17392A5  = 0.02936 
0.62110A1-1.59420A3+0.52294A5  = 0.02731 
0.17391Ax+0.52293As-1.6S986A5  = -0.00799 
These  two  sets  of  equations  may  be  solved  for  the  A's,  to  give 


Aq  = 0.057659  and  A^  = 0.022833 
A2  =-0.031776  A3  =-0.006525 
A4  = 0.002552  A5  = 0.005213 
A6  = 0.041474 


Likewise,  the  values  for  the  Am+jc  coefficients  are  obtainable  from  the 


relations 

A7  = 


4 


// 


/ 

^7,7 

(&, 

-^7 -A 

> A 

X ~ A + 

'A*  A‘ 

/ 

y^9j9 

A 

“■^*9.0  ^ 0 

"^92  Az 

\ 

\) 

/ 

^n,n 

a 

A> ' 

''^6 

/ 

U, 

-'“s.A  ■ 

■'«VA') 

Ale-  y.  (t'/o  V*/<y  A s 

^ ^ ~ (^iz~ ^ ^s) 


so  that  Ay  = 0.187536  and  Ag  = 0.110043 

Ag  = 0.107867  A10=  0.079285 

An=  0.056964  Ai2=  0.065819  . 

It  is  worthwhile  to  make  the  following  observation  about  these  results 
at  this  juncture.  Let  these  just-determined  values  of  An  and  Aju+fc  be 
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inserted  into  Ea.  (57)  of  Part  I in  order  to  find  the  corresponding  value  of 
S**.  Now  if  the  so-determined  S**  function  is  \xsed  in  Eq.(48)  of  Part  I,  then, 
because  of  the  physical  significance  attached  to  the  e quantities  that  was 
pointed  out  in  the  theoretical  development,  it  follows  that  the  exact  values 
of  the  drag  coefficient  will  be  obtained  which  pertains  to  a duct,  the 
radial-velocity-component  parameter  of  which  is  defined  by 

rp,)(i^e4-e) 

and  which  is  represented  by  a trigonometric  series  in  cos(nG)  and  sin  (nG), 
with  n taking  on  integral  values  from  0 to  six  in  the  particular  example  now 
under  examination. 


DETERMINATION  OF  THE  BEST  DUCT  CONTOUR  UNDER  CONDITION  (a) 
TOGETHER  WITH  CONDITION  (b)  OR  (c) 


6.  Calculation  of  the  Coefficients  Hn  and  Kjj 


The  values  of  the  coefficients  Hn  and  can  be  determined  by  use  of 
Ecs.  (65)  of  Part  I,  once  the  coefficients  and  Cn  appearing  therein  have 
been  evaluated.  This  evaluation  for  the  Bn  and  Cn  coefficients  is  carried 
out  by  having  recourse  to  Ecs.  (64)  of  Part  I.  Making  use  of  the  specific 
date  valid  for  this  example,  it  is  found  that 
B0  = 3.14159  Al  + 4.71239  /^2 
B1  = -1.90984 
B2  = 0 

Bj  = -0.21220  A2 


B4  = 0 

B5  = -0.07639  A 2 
B6  = 0 
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while 


C0  = 0 


C4  = -0.75  A 2 


Ci  = 2\1  + 5/^2 


C5  = 0.4/^  + 0.6^2 


C2  = -1.5/1  2 


Cg  = -0.5 2 


C3  = 0.66667  A 1 + A2 

In  addition  to  these  quantities  required  for  insertion  into  Eq.  (65)  of 
Part  I)  it  is  also  necessary  to  know  the  values  of  c and  o'  appearing  therein. 
The  values  of  c and  c'(  may  be  obtained  by  use  of  Eqs.  (66)  of  Part  I,  and  the 
pertinent  calculations  result  in  the  following  values: 


c0  = 1 - 1.5  A0  - 0.5625  Aq  - 0.88108 

cx  = 1 — 0.75A1  - 0.28125A-L  = 0.97645 

c2  = 1 — 1.0 3125 A2  = 1.03277 

c3  = 1 - 1.03125Aj  = 1.00673 

c4  = 1 - 1.03125A4  = 0.99737 

c5  = 1 - 1.03125A5  = 0.99462 

c6  = 1 - 1.03125A6  = 0.95723 

while 


C1 


c3  = 


c4 


C6 


1 

1 

1 

1 

1 

1 


1.03125A7  = 0.80660 
1.03125Ag  = 0.88652 
1.03125Ag  = 0.83876 
1.05125Aio=  0.91824 
1.03125AU=  0.94126 
1.03125A12=  0.93212 


Finally,  the  coefficients  br,i  and  bi,r>  which  are  defined  as  part  of 


Eqs.  (66')  in  Part  I,  may  be  seen  to  have  the  values: 
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0 


*0,1  = 2 
bo, 2 = 1 

bQ  3 = 0.66687 

bo, 4 = °»5 
bo, 5 = 0.4 
bQ)6  = 0.53335 
and 

b4,0  = 0 
b4,i  = -0.13333 
b4,2  = -0.33353 
b4>3  = -0.85714 

b4,5  = 1.11111 

b4,6  = °*6 
and  likewise 

bl,0  = 0 
bl,2  = 1.33333 

bl,3  = °*75 
bl,4  = 0.53333 

b^s  = 0.41667 
bl,6  = 0.34286 
and 

b5,0  = 0 
bs^].  = -0.08333 
b5,2  = -0.19048 
b5,3  = -0.375 
b5>4  = -0.88889 
b5^6  = 1.09091 


b2,0 

b2,l 

b2,3 

b2,4 

b2,5 

b2,6 

b6,0 

b6,l 

b6,2 

b6,3 

b6,4 

b6,5 

b3,0 

b3,l 

b3,2 

b3,4 

b3,5 

b3,6 


-0.66667 

1.2 

0.66667 

0.47619 

0.375 

0 

-0.05714 

-0.125 

-0.22222 

-0.4 

-0.90909 

0 

-0.25 

-0.8 

1.14286 

0.625 

0.44444 
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In  order  to  get  a first  approximation  to  the  sought  H and  K coefficients, 
one  may  now  make  use  of  the  above-determined  quantities,  and  guess,  as  a first 
set  of  values  for  L^,  L^',  Lr  and  Lr',  to  be  inserted  into  Eos*  (65)  of  Part  I, 
the  common  value  of  zero.  The  equations  for  the  determination  of  the  first 


approximate  values  for  the  H's  then  turn  out  to  be  simply: 

0. 83103  H0  = 

Bo+0- 47746  ^ 0,46500  [ci+0,47746(0-08295Ho-0,07336H2-0o01057H4-0. 00376^3 
+0  o 06091  [c  3+0  „ 47746  ( 0 c 027  65Ho+0 . 13205H?.-0 . 06796H4-0 . 01463H6J 
+0  „ 02421  ]05+0  o 47746  ( 0 „01659Ho+0 . 0524GH2+0 .0880914-0  ,05983145] j 
1,03277H2  - 

B2+O. 47746  |^-0o  15500  [ci+0„03960Ho-0„03503H2-0  .0050514-0. 00179He] 

+0 . 14564^03+0 ,01320H0+0o06505H2-0o03245H4-0 „00698H6] 

+0  o 02882  [C5+0 . 00792Ho+0  ,02502H2+0 -04206H4-0  .02357H^jj 
0., 99737  H4  = 

B4+0 ,47746* -0  ,03100 [ci+0.03960Ho-0 ,03503H2-0 „00505H4-0 -00179H6] 

-0 . 10403  |c 3+0 , 01320Ho+0 „06305H2-0 „03245H4-0  .00698H6] 

+0 . 06724  [c 5+0 , 007 92Ho+0  „02502H2+0  .0420614-0  <,02357H6]j 
0 ,95723  H6  =■ 

Bg+0  o 47746  j"-0  ,01328  [ci+0 ,03960Ho-0 .03503H2-0  ,0050514-0  „00179H^[ 

-0  a 02697  [C3+O „01320Ho+0  o06305H2-0 . 03245H4-0 .00698H^[ 

-0 . 05502  [C5+O  „ 00792Ho+0  .02502H2+0  „04206H4-0  „0-2357Hef^ 

0.97645  Hx  = 

B 2+0-47746  j 0ol6550  £2+0  , 47746(0, 25005fli-0.08629H5-0.01085!I5j] 

+0  „ 04605  [c4+0  -47746(0, 10002H1+0 , 12528fls-0  -05063^)] 

+C  ,02421  [cg+0,47746(0-06430Hi+0-04794H3+0-06214H5)j'l 
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1.00673  H5  = 

B3+0. 47746  j -0.09930 [c2+0.11939Hi-0.04120H3-0-00518H5] 
+0.09868  |04+0 .04775Hi+0 .05886H3-0.02417H5] 
+0 .03138  [06+0 .03070Hi+0 . 02289H3+0  .02967HsJJ 


0.99462  Hg  = 


B5+0. 47746  |-0.02364[C2+0. 11939^-0. 04120H3-0.00518H5] 
-0.07675[c4+0.04775Hi+0.05886H3-0.02417H5] 
+0. 07703  [c6+0 .03070Hi+0 .02289H3+0  .02967H5]j 


These  equations  may  be  reduced  in  complexity,  therefore,  to  give  just  the 
following 

0 . 87169Ho+0 „00505H2+0 .OOI88H4+0 .00100H6=BQ+0 .2220201+0 .0386303+0 .0115605 
0 .00190Ho+l .02546H2+0 .00131H4+0 .00074H6=B2-0  .0740101+0 .069540 3+O .01376CS' 

0 .00100Ho+0 .OOI8IH2+O .99434H4+0 .00054H6=B4-0 .01480Ci-0 .0496703+0 .05210Cs 
0.00063Ho+0.00125H2+0.00065H4+0.95638H6=B6-0.00634Ci-0.01288C3-0 .0262705 
and  likewise 

0 . 965 62H1+0 . 00171H 3+0 . 000 6OH5  =Bl+0 .0790202+0 .02194C4+0 .0115606 

0 . 00295H1+1 .00167H3+0 .00045Hs=B3-0 .0474102+0 .04712C4+0 .01498C6 

0 .00197H1+0 .00086H3+0 . 99258H5=B5-0 .0112902-0 .0366404+0 .0367806 . 

Before  application  of  the  pertinent  boundary  conditions,  the  value  of 

A j and  A 2 are  unspecified  parameters,  so  that,  at  this  stage,  the  solution 

for  Hn  and  will  be  of  the  form  Hn  = Hn* A 1 + 2 anci  Ai+Kfc*  A 2* 

Thus,  in  proceeding  to  the  next  step  in  finding  the  first  approximate  value 

(1)  _ *(1)  . **(1)  . 

for  Hn  one  may  look  for  it  in  the  form  Hn  - nn  A 1 + Hn  A 2>  vhere 
the  superscript  "one"  denotes  that  the  values  to  be  obtained  are  but  the  first 
approximations.  Upon  separation  of  the  reduced  equations  just  derived  above 
into  the  part  which  is  multiplied  by  A 1 and  the  part  which  is  multiplied  by 
?[2t  therefore,  it  will  be  found  that  the  terras  linked  to  A 1 give 
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0.87169Ho  +0.00505H2  +0.00188H4  ^ ^ +0»00100Hg  3.61601 

0.00l90Ho*(l)  +1-02546H2*(1)  +0.00131H4*(1)  +0.00074^*^  = -0.09616 

0„00100Ho*(1)  +0o00181H2*{;,)  +0.99434H4*(l)  +0.00054H6*(l)  =-0.04987 

0.00063H0*^  +0.00125H2*^  +0.00065H4  +0.95638H6*^1^  = -0.03177 

#(  1 \ 

and  solution  of  this  simultaneous  set  of  four  eouations  gives  the  H v ' values 


as 

*(1) 

H0  = 4.14828 

*(1) 

H2  = -0.10146 

- -0.05414 

*(l) 

Hg  = -0.03578 

while,  because  Xl  does  not  appear  in  the  expressions  for  Bi,B3,B5,C2>C4,  and 


Cg,  if  follows  that 


Hi 


*(1)  = H_*(l)  _ Ur*(l)  = 


= h5 


= 0. 


In  similar  manner  it  will  be  found  that  the  terms  multiplied  by  \ 2 must 


obey  the  relations: 

**(1)  **(1)  **(1)  **(1) 

0.87169Ho  +0.00505H2  +0.00168H4  +0.00100H6  = 5.42402 

**(1)  **M)  *#(l)  **(l) 

0.00190Ho  v y+1.02546H2  V '+0.00131H4  v +0.00074H6  = -0.14423 

1 ) **(1)  **(l)  **(l) 

0.00100H0  'J‘y+0o00181H2  v +0.99434H4  v '+0.00054Hg  = -0.07481 

o.oooesHo^^Woo^H^^WoooesH^^+o.gsessHe**^  = -0.04766 


and  likewise 

0,96562H^  Hf^+0„00171H5**^+0.00060H5  ^ = -2.05064 

0.00295H1**(1^+1.00167H3'!W^140.00045H5^^1^  = -0.18392 
0.00197H1^(1)+0.00086H3^(1)-K).99258H5^(1)  = -0.05036 


and  upon  solution  of  these  two  sets  cf  simultaneous  equations,  it  is  found  that 
**(1)  = 


H**U)  . 6.22242 
no 


h2 


-0.15218 


**(l) 

H4  v 7 = -0.08122 


H 


■**(1) 


= -0.05368 
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and 


= -2 o 12365 


Hs**^  = -0.17736 


HS 


**(D  = -0, 


04637 


By  use  of  the  approximation  that  Lj_  = Lj_'  = Lr  =Lr  =0  still,  it  is  likewise 
seen,  by  use  of  Eos.  (65)  of  Part  I,  that  the  first  approximate  values  for  the 
K's  are  determined,  in  the  present  instance,  by  use  of  the  following  simple  ex- 
pressions, which  have  actually  been  employed  already  in  the  determination  of 


the  H-values: 

Kx  = 1.23977  [c^O. 47746(0. 08295Ho-0.07336H2-0.01057H4-0.00376H6^ 

K*  = 1.12516  [C3+O. 47746(0. 02765Ho+0.13205H2-0.06796H4-O„01463H6j 
K5  = 1.06241  [C5+O. 47746(0. 01659H0+0.05240H2+0.08809H4-0.05983H6)| 

K2  = 1.12801  [c2+0 .47746(0.2 5005H]_-0 . 0 8 62 9H  3-0 . 010 8 5H^)j 
K4  = 1.08904  [C4+O. 47746(0. 10002Hi+0.12528H3-0.05063H5j 
K6  = 1.07282  [Cg+0. 47746(0. 06430H1+0.04794H3+0.06214H5^ 

As  was  done  in  the  case  of  the  H's,  the  K's  are  to  be  separated  into  two 
parts  linked  to  the  A ^ and  A 2 parameters,  and  if  the  superscript  "one"  is 
again  used  to  denote  the  first  aporoximation  values,  the  solutions  of  the 


above  equations  are  found  to  be; 

*(1) 

Kx  = 2.68804 

K^1)  = 0.80677 


* 


Kg  (1)  =0.45532 

*(1)  *(1)  *(1) 

while  K2  = K4  = K6  = 0 


and,  likewise,  the  solutions  for  the  part  of  the  K's  which  is  linked  to  the 


parameter  are  found  to  be,  for  the  first  approximation: 

**(1) 

Kx  = 4.03207 


**(1) 


k5 


-JHf- 


*(D 


= 1.21015 
= 0.68373 
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**(1) 

% = -1.96949 

K4^(i)  = _o.  93733 

%**^  = -0.61219  . 

Now  that  the  first  approximation  values  for  the  H's  and  K's  are  at  hand, 
it  is  possible  to  compute  the  corresponding  first  approximation  values  for  L 
and  L',  according  to  Eqs.  (66')  of  Part  I.  The  general  equations  for  deter- 
mining the  successive  L's  and  L^s  are  thus 
L0  = -0.04567  Hi  +0.00145  H3  -0.00042  H5 

-1. 57080(-0 . 37507K1+0 . 11004K2-0 .07191K3+0 .03964K4-0 .02279%+0 .02194%) 

% = -0.11532Ho+0.03531H2-0.00039H4-O.00251H6 

-1.57080(-0.14672K2+0.03090K3-0.04228K4+0.02373%-0.02257K6+0. 13754%) 
L2  = -0.02537H!+0.00679H3-0.00069H5 

-1 . 57080(0 .12502K1-0 .12944K3+0.05236K4-0 .02713%+0 .02468%+0 .05502%) 

L3  = -0.01281Ho+0.03305H2-0.00260H4-0.00512H6 

-1 .57080(-0 .04688K1+0 .08803K2-0 .0906114+0 .0687a%-0 .02925%+0 .03596%) 
L4  = -0.00345H1+0.00666H3-0.00528H5 

-1. 57080(0. 025C0%-0.03668K2+0.09246K3-O.06329%+0. 0394 9%+0. 01932%) 

L5  = -0.00461Ho+0.00413H2-0.00258H4-O.04132H6 

-1 . 57080  ( -0 .01563%+0 .02096K2-0 .04045K3+0 .07048%-0 .07130%+0 .01316%) 

L6  = -0 .00138Hi+0 . 00081H3-0 .00526% 

-1.57080(0 .01072%-0 .01375K2+0 .02397%-0 .03171K4+0 .05178%+0 .01097%) 
and  the  similar  array  for  the  primed  L*s  has  the  specific  form 

= -1 .57080(-0 .11532Ho-0 .02118H2+0 .00163H3+0 .00034H4-0 .00043H5+0 .00237H6 
+0 .02283Hi) -0 . 12227K2-0 .01198%-0 .00398% 

I 

L2  = ^1 .57080 (0 .05766Ho-0 .03044HX-0 .00522H3-0 .OOO85H4+O .00099Hs-0 .00518% 

-0 .01589H2)  -0  .20837%-0  .11218%-0 . 00  749% 
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Lj  = -1  .57080(-0.03844Ho+0 . 01712Hi+0  .03813H2+0  .0021914-0  .00195H5+0  .00922% 

-0  .00217%) -0 .11444*2-0 .08090*4-0 .00813K6 
L4  = -1 .57080  (0 .02883Ho-0  .01213^-0  .0211814+0  .OO746H3+0 .00463H5-0  .01659Hg 
+0 .00064H4) -0 .02834K1-0 .11007*3-0 .05767*5 
L5  = -1 . 57080  ( -0 .02306Ho+0 .00951Hi+0 .01513H2-0  .00408H3-0  .00284H4+0  .03770H6 
+0 .00104H5) -0 .01447*2-0 .08026*4-0 .06636% 

Lg  = -1.57080(0. 01922H0-0.00783H1-0.01192H2+0.00290H3+0.00153H4-0.00569H5 
+0  .00691%)  -0 .01133*1-0 .01332*3-0 .05743% . 

Treating  these  L and  L'  values  in  the  same  way  that  the  H's  and  K's  were 
handled,  by  splitting  them  up  into  two  parts,  the  evaluations  will  be  sought  in 
the  form 

L = L*/\x  + L**A2  and  L'  = L'*^i  + L'^2- 
The  first  approximations  may  now  be  obtained  directly,  therefore,  by  sub- 
stituting the  above-determined  first  approximation  values  for  H*,H**,K*,  and  *** 
into  the  equations  for  L and  L'  just  set  down  above.  The  numerical  results  ob- 
tained when  this  procedure  is  carried  through  are  found  to  be 


L*  = 1.69114 

"0 

= 

3.05334 

% = -1.39323 

T** 

L1 

= 

-2.62770 

= -0.34441 

^2 

ss 

-0.19215 

L5  = 0.04695 

L3 

= 

0.18123 

* 

** 

L4  = -0.17742 

L4 

= 

-0.30606 

** 

L5  = 0.08992 

L 5 

— 

0.23447 

Lg  = -0.11272 

T** 

Le 

= 

-0.24472 

and  likewise  it  is  found  that 
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1.45339 


Li  = 0.74822  Ia  = 


1* 

1*2 

= 

-1.03263 

L 2 

- 

-1.65188 

t* 

'#* 

L3 

= 

0.25727 

L3 

= 

0.74846 

i* 

t** 

^4 

= 

-0.38338 

l4 

= 

-0.61327 

• ■JHf 

0.40684 

4 

= 

0.15455 

l5 

— 

i* 

t-JHt- 

l6 

-0.19400 

L6 

= 

-0.31672 

If  the  method  of  successive  approximations  is  now  carried  one  step  farther 
by  use  of  the  just  found  values  for  L and  L’,  it  will  be  seen  that  the  two  sets 
of  simultaneous  equations,  from  which  the  second  approximation  values  for  H and 
K may  be  determined,  become  now 

0.87169Ho*(2)  +0.00505H2*(2)  +0.00188H4*(2)  +0.00100H6*(2)  = 3.82903 

0.00190Ho*(2)  +1.02546H2*(2)  +0.00131H4*(2)  +0.00074Hg*(2)  = -0.13944 
0.00100Ho*(2)  +0.00181H2*(2)  +0.99434H4*(2)  +0.00054H6*(2)  = -0.07224 
0.00063Ho*^2^  +0.00125H2*^2^  +0.00065H4*^  +0.95638Hg*(2)  = -0.04600 
and  for  the  set  with  odd  subscripts  one  has 
0.96562H!*^  +0.00171H3*^  +0.00060H5*(2)  = -0.16932 

0.00295%*^  +1.00167H3*^  -K).00045H5*^  = 0.00354 

0.00197Hi*^  +0.00086%*^  +0 .99258H5*^  = 0.01237  . 

Upon  solving  these  sets  of  equations  it  is  found  that  the  second  approxi- 
mation values  for  the  H's,  at  least  that  part  of  the  H which  is  related  to  the 

Xi  parameter,  are  as  follows; 

H0*(2)  = 4.39371  = -0.17537 

h2*^  = -0.14398  H3*(2)  = 0.00904 

H4*(2)  = -0.07678  H5*(2)  = 0.01280  . 

H6  w = -0.05075 

By  following  the  same  procedure  as  utilized  in  the  previous  approximations,  the 
corresponding  values  of  K may  be  obtained  from  use  of  these  H values  just  set 


down  above,  and  upon  carrying  out  the  pertinent  calculations  it  is  found  that 
the  second  approximations  for  the  K's,  at  least  the  part  linked  to  the  Ax 
parameter,  turn  out  to  be 


*(2) 

Kl 

= 2.80785 

#(2) 

K4  K ' = -0.05647 

K *(2) 
k2 

= -0.15638 

Kg*^  = 0.47491 

*(2) 

= 0.84134 

*(2) 

Kg  = -0.02887  . 

In  an  analogous  fashion  the  H values  which  are  linked  to  the  A 2 parameter 

are  determined  by  solution  of  the  array  of  equations 

**(2)  *#(2)  **(2)  #*(2) 

0.87169Ho  +0 .005051*2  +0.00138K4  +0.00100Hg  = 5.81265 

0 .00190Ho**  ^ 2 ^ +1 .02546H2**  ^ 2 ^ K> . 00131H4**  ^ 2 Wo0074Hg^  2 ^ = -0 . 171 82 

0.00100Ho^2)+0.00131H2^2^+0.99434H4^2^+0.00054H6**^  = -0.11489 

0 .00063Ho**(2)  +0 .00125H2** ( 2)  +0  .00065H4** ( 2)  +0 . 95638H0** '2) = -0 .07892 


and 

-**(2)  **(2)  **(2) 

0.96562%  V '+0.00171H3  +O.OOO6OH5  = -2.36698 

-**(2)  #*(2)  **(2) 

0.00295H1  ;+1.00167H3  +0.00045%  = -0.15817 


0.00197H- 


**(2) 


**(2) 


+0.00086H3  +0.99258% 


*#(2) 


= -0.02028. 


The  solution  of  these  equations  will  be  found  to  give  the  following  values 


##(2) 

for  the  H values: 


H, 


**(2) 


= 6.66966 


H2**(2)  = _o. 17969 

H4**^  _ _q#12108 

H6^(2)  _ ^^Qgggo 


**(2) 

% = -2.45098 

**(2) 

H3  = -0.15068 

**(2) 

H5  K ' = -0.01544 


By  reuse  of  these  values  in  the  way  that  it  was  done  for  the  first  approxi- 
mate) 

mation,  it  will  be  found  that  the  corresponding  evaluations  for  the  K 


values  are 
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Ki**(2)  _ 4.26093 
K2**(2)  = -2.22738 
K3**(2)  = 1.31258 
K4^(2)  _ _x. 02963 

Ks**(2)  = 0.73521 
Kg** (2)  = -0.66007 


There  appears  to  be  little  reason  to  proceed  to  any  further  degree  of 
approximation  by  continuing  the  above  process,  so  that  consequently  it  may  now 
be  taken  for  granted  that  the  H and  K values  pertinent  to  the  case  in  hand,  have 
the  following  evaluations: 


while 


Ho 

H1 

h2 

h3 

h4 

h5 

h6 


4.39371  + 6.66966  /l x 

-0.17537  Xl  “ 2.45098  A2 
-0.14398  Ax  - 0.17969  \2 
0.00904  Ax  - 0.15068  A 2 
-0.07678  Ax  - 0.12188  A2 
0.01280  Al  - 0.01544  A 2 
-0.05075  Al  - 0.08660  >2 


*1 

*2 

*3 

K4 

K5 

*6 


2.80785  Ai  + 4.26093  A2 

-0.15688  \ - 2.22738  A2 

0.84134  Ai  + 1.31258  A2 

-0.05647  Ax  - 1.02963  A2 

0.47491  Ai  + 0.73521  A2 

-0.02387  Al  - 0.66007  A 2 


The  functional  relationship  which  describes  how  the  slope  of  the  external 
contour  of  the  duct  must  behave  under  the  thus-far  imposed  general  conditions 
may  now  be  written  down  by  aid  of  the  A,  H,  and  K values  that  have  been 
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determined  above.  If  the  sought  fj  -function  is  broken  up  into  the  part  linked 
to  the  A]i -parameter  and  the  part  linked  to  the  X 2-parameter  by  use  of  the 
definition  that 

f})  = 

then  it  follows,  from  use  of  Eqs.  (56)  and  the  one  preceding  it  in  Part  I,  that 
the  respective  parts  of  the  fj -function  have  the  fora: 

rj^z)  = J^25z  + 0.12096-0.00191  cos  9 + 0.00218  cos  20 

-0.00005  cos  50  -0.00010  cos  40  +0.00005  cos  59  -0.00100  cos  6© 

+0.25142  sin  0 -0.00824  sin  20  +0.04555  sin  50  -0.00214  »in  40 

+0.01292  sin  50  -0.00091  sin  60 

and 

= ;--z  -b  + 0.18562  -0.02672  cos  0 + 0.00275  cos  20  +0.00047  cos  30 
12'  1+0.252' 

-0.00015  cos  49  -0.00004  cos  50  -0.00171  cos  60  +0.38155  sin  0 

-0.11703  sin  20  +0.06760  sin  30  -0.03898  sin  40  +0.02000  sin  50 

-0.02074  sin  60. 

7.  Determination  of  the  Ai  and  A 2 Parameters 

The  stage  has  now  been  reached  in  the  development  of  the  computations  where 
it  is  possible  to  apply  the  final  specific  perimetral  conditions  which  were 
selected  for  study  in  the  Introduction.  These  cases  will  now  be  restated, 
and  the  final  step  in  the  determination  of  the  duct's  slope  will  then  be 
carried  out. 

Case  (a) : In  this  case  the  sole  condition  to  be  imposed  is  that  the  duct 

should  close,  i.e.,  it  is  merely  assumed  that  the  "end-points"  and  N2 
at  the  lip  and  at  the  tail  end  of  the  sought  duct's  meridional  contour  are 
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to  belong  to  the  frustum  of  a cone  section  which  constitutes  the  inner  duct  whose 

Ro  r 

large  end  and  small  end  have  radii  in  the  ratio  -=  = 1 + 8 * r . where  Ri  is  here 

Rx  Rx  A 

equal  to  the  radial  coordinate  of  the  Nx  lip-point  and  R2  is  equal  to  the  radial 
coordinate  of  the  N2  tail-point  of  the  sought  external  contour. 

In  this  instance)  therefore)  the  appropriate  conditions  on  the  constants 
Al  and  X 2 are  that  ^2  i®  to  be  taken  equal  to  zero,  while  x is  found  by 
making  use  of  Eq.  (68)  of  Part  I,  which,  in  the  present  notation)  beoomes 


l , EliEL 

6 R, 


" b r ^ ' ^if  ^ y 


Upon  use  of  the  appropriate  concrete  data  applying  to  the  case  in  hand)  it 
is  found  that 

J \ = &1H05 

and  therefore  the  value  of  A 1 turns  out  to  be  X r --?✓  ■■ 

* ' / 3.  //+« 05 

In  the  Normalizing  Mach  Number  Case,  therefore)  where  one  has  that  B = 1 
and  that  0 = 0.25,  the  corresponding  value  of  A 1 i®  ^ 1 = 0.24094.  The 
slope  function  which  applies  in  this  specific  case  then  becomes 

1 (r)  - x»0'"25z  +0*02915-0-00046  cos  9 +0.0052  cos  2©  -0.00001  cos  5© 
-0.00002  cos  4©  +0.00001  cos  5©  -0.00024  cos  6©  +0.06055  sin  © 

-0.00198  sin  29  +0.01044  sin  5©  -0.00051  sin  4©  +0.00511  sin  50 

-0.00022  sin  6©  . 

The  actual  R coordinates  of  the  sought  duct  are  then  given  by  the  integral, 
which  takes  the  form 
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j n(s)ds  = = 0.96336  loge  (1+0.25^  +0.02913y  -0.00044  sin(- £%■) 

° + 0.00025  sin  (—^)  -0.00000  sin  {Vf)  -0.00000  sin  (~*j 

+ 0.00000  sin  00004  sin  (2 77^) 

+ 0.05782  ^1— cos  (ZlS^  -0.00094  [l-cos 

+ 0.00332  1-cos  (Tfy^  -0.00012  1-cos  ('ApC^ 

+ 0.00059  1-cos  (^p^J]  -0.00003  1-cos  (2/^jj  . 

It  will  be  instructive  to  find  out  how  the  two-dimensional  solution  compares 
with  this  result.  In  other  words  a computation  will  now  be  made,  using  the  same 
perimetral  condition  as  employed  in  the  above  example,  but  this  time  the  simple 
Ackeret  formula  will  be  utilized  for  relating  the  local  pressures  to  the  local 
slopes  of  the  sought  duct's  contour.  Under  these  circumstances  the  local  slopes 
of  the  contour  are  described  by  the  formula 

n = Ai. 

I 1+0 .25 z 

and  consequently  the  value  of  Ax  found  to  be  Ax  = 0.33505,  because  here 


it  is  true  that 


j = 2*23348 

Ja 


The  formula  for  the  "two-dimensional”  duct  coordinates  thus  becomes 


r-rx  n 

BRx  / 


n 

■ fj(s)ds  = 1.34020  loge  (1+0.25^). 

SO 


R— Rl 

The  radial  coordinates,  — , which  apply  to  the  duct  in  question,  which 

has  minimum  drag,  with  no  perimetral  condition  imposed  other  than  the  one  of 
closure,  and  in  the  Normalizing  Mach  No.  Case,  have  been  computed  for  certain 
selected  abscissa  stations,  y , along  the  duct's  length.  The  results  for  both 
the  three-dimensional  and  the  two-dimensional  cases  are  listed  in  Table  XI 
below,  while  the  plot  of  the  meridional  line  now  obtained  is  shown  in  Fig.  5 
for  both  cases  also. 
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TABLE  XI 


DUCT  COORDINATES  FOR  CASE  OF  NO  PERIMETRAL  CONDITION.  SAVE  THAT 
OF  CLOSURE  - NORMALIZING  MACH  NO. 


y 

K- 

BR, 

.RzRi 

B/?, 

By  Use  of  Formula  (11) 

By  Use  of  Formula  (12) 

0 

0 

0 

0.4 

0.11113 

0.12773 

0.8 

0.22363 

0.24434 

1.2 

0.33155 

0.35161 

1.6 

0.43540 

0.45094 

1.8 

0.48555 

0.49796 

2.2 

0.58181 

0.53734 

2.4 

0.62814 

0.62989 

2.8 

0.71372 

0.71115 

3.0 

0.75000 

0.75000 

In  the  Small  Conical  Diffusor  Angle  Case  the  value  of  will  be 

changed  from  what  it  was  above,  but  the  same  procedure  will  otherwise  be 
valid.  In  this  new  case,  therefore,  one  has  that  B = 1.66667  and  p = 0.15; 
so  that  one  may  now  proceed  as  in  the  former  analogous  case  to  obtain 


0.45 

3.11405 


0.14451 
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DUCT  CONTOUR  FOR  LEAST  DRAG  IN  CASE  OF  NO  PERIMETRAL 
CONDITION  SAVE  THAT  OF  CLOSURE  - NORMALIZING  MACH  NO.  CASE 


DUCT  CONTOUR  FOR  LEAST  DRAG  IN  CASE  OF  NO  PERIMETRAL 
CONDITION  SAVE  THAT  OF  CLOSURE  - small  conical  DIFFUSOR  ANGLE 

CASE 


LEGEND 

MERIDIONAL  CONTOUR  OBTAINED  BY  ASSUMPTION  THAT 

PRESSURE  COEFFICIENTS  ARE  RELATED  TO  DUCT 

GEOMETRY  BY  MEANS  OF  ~EQ.  (47)  OF  PART  J 

MERIDIONAL  CONTOUR  obtained  BY  ASSUMPTION  THAT 

PRESSURE  COEFFICIENTS  Jjffe  ARE  RELATED  70  DUCT 
PLANAR  ELEMENTS  THROUGH ~ ACNERET'S  TWO' DIMENSION' 
AL  FORMULA 
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nr 


and  thus  the  corresponding  description  of  how  the  duct  contour  behaves  is  provided 
by  the  function 


= 0.57804  loge  (1+0*25^)  +0.017487  -0.00027  sin  (—*) 


+ 0.00015  sin  -0.00000  sin  ( tty)  -0.00000  sin  (^*f) 

+ 0.00000  sin  -0.00002  sin  (2 it}) 


+ 0.03469 

+ 0.00199 
+ 0.00036 


H (¥)] 


1-COS 

1-cos 


-0.00057 

-0.00007 

-0.00002 


1-cos 

1-cos 

1-cos 


(¥*)] 


(ii') 


By  use  of  the  Ackeret  formula  under  the  same  conditions  one  would  obtain  on 
the  other  hand  that 


X 


1 = 


0.45 

2.23848 


0.20103 


and  thus  the  ordinates  in  this  case  would  be  given  by  the  expression 

R_R, 

= 0.80412  loge  (1+0 .25 z) . (12') 

R— Rl 

The  radial  coordinates,  ■ ■,  which  apply  to  the  duct  in  question,  which 

BRl 

has  minimum  drag,  with  no  perimetral  condition  imposed  other  than  the  one  of 
closure,  in  the  Small  Conical  Diffusor  Angle  Case,  have  been  computed  for  certain 
selected  abscissa  stations,  z,  along  the  duct's  length.  The  results  for  both 
the  three-dimensional  and  the  two-dimensional  cases  are  listed  in  Table  XII 
below,  while  the  plots  of  the  meridional  line  obtained  by  use  of  these  data  are 
shown  in  Fig.  6 for  both  instances. 
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TABLE  XII 


DUCT  COORDINATES  FOR  CASE  OF  NO  PERIMETRAL  CONDITION. 
SAVE  THAT  OF  CLOSURE  - SMALL  CONICAL  DIFFPSOR  ANGLE 


y 

/?-/?, 

BR, 

R-f?, 

BR, 

By  Use  of  Formula  (11') 

By  Use  of  Formula  (12') 

0 

0 

0 

0.4 

0.06669 

0.07664 

0.8 

0.13418 

0.14661 

1.2 

0.19892 

0.21097 

1.6 

0.26123 

0.27056 

1.8 

0.29133 

0.29878 

2.2 

0.34909 

0.35241 

2.4 

0.37689 

0.37794 

2.8 

0.42825 

0.42669 

3.0 

0.45000 

0.45000 
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Caae  (b):  "No  Increase  in  Area  Over  that  Occupied  by  Frustum-of-Cone  Basic 

Shape".  In  addition  to  the  closure  condition  employed  in  Case  (a),  it  is  now 
further  postulated  that  the  area  enclosed  between  the  sought  contour  and  the 
x-axis  is  to  be  equal  to  the  area  which  would  be  enclosed  between  this  x-axis 
and  the  straight  line  segment  joining  the  points  N^  and  N2» 

Under  this  general  perimetral  condition,  the  Al  and  A 2 parameters  will 
have  to  be  such  as  to  permit  Eqs.  (68)  and  (69)  of  Part  I to  be  satisfied.  In 
the  present  instance,  therefore,  it  will  be  required  that  the  specific  form 
of  Eqs.  (68)  to  be  satisfied  should  be 


q/t) 

0 

while,  because  the  constant  C appearing  in  Eq.  (69)  takes  on  the  specific 
B L2 

value  C = £ . — “r  when  the  area  enclosed  by  the  sought  contour  is  not  in- 

2 B2R-J2 

creased  over  the  area  enclosed  by  the  straight  line  chord  between  the  end- 
points Ni  and  N2  and  the  x-axis,  it  follows  that  in  this  present  instance 
Eq.  (69)  reduces  to  just 


fa 


X , dy  + ‘dy  ( - -4  * • 

•sq  6 ^ 

In  the  Normalizing  Mach  No.  Case  for  which  B = 1,  and  (J  = 0.25,  therefore, 
the  evaluated  forms  of  the  above  perimetral  conditions  become 


rl, 


5.11405  Al  + 4.37629  A2  = 0.75 
4.96539  Ai  + 5.06409 A2  = 1*125 

from  which  it  follows  that  the  pertinent  values  of  Al  and  A 2 ar® 

Ai  = 0.18879  and  A 2 = 0.05704  . 

Likewise,  in  the  Small  Conical  Diffusor  Angle  Case,  for  which  B = 1.66667 
and  6 = 0.15,  it  is  readily  seen  that  the  evaluated  forms  of  the  above-stated 
perimetral  conditions  are 
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3.11405  Ax  + 4.37629  A2  = 0.45 
4.96539  Ax  + 5.06409  X 2 = 0.675 
from  which  it  follows  that 

Al  = 0.11328  and  A2  = 0.02222. 

Upon  insertion  of  the  \ -values  into  the  expression  for  the  sought  slope 
function  it  is  then  found  that 

In  the  Normalizing  Mach  No.  Case: 

1 <»>  = liSif  + MS*  -0-00155  cos  (**-) 

+0.00051  cos  .00001  cos  (vy  )-0.00004cosppy) 

+0.00000  cos  ^^|i)-0.00025  cos(2^^ )+0.06160sin 
-0.00589  sin  (^p)+ 0.01068  sin(  ) -0.00185 sin 

+0.00318  sin  (^)-0 .00094  sin (2;^)  . 

The  actual  R coordinates  of  the  sought  duct  may  then  be  obtained  through 
integration,  with  the  result  that 

= 0.16252  logg  (1+0.25^  +0.17780f -0.00129  sin  (£Jl) 

+0.00024  sin  0.00000  sin  {fry)  -0.00001  sin  (~g~) 

+0.00000  sin  (^p^-0 .00004  sin  (2 Ty)  +0.05882  [l-  cos  (^~) ] 

-0.00281  [l— cosfnpj  +0.00340  [l-cos  (^) ] -0 .00044 |l-co s 

+0.00061  [l-cos(^j] -0.00015  [l-cos  (2^  (13) 

In  the  Small  Conical  Diffusor  Angle  Case: 

The  analogous  formula  for  the  slope-function  is 
fj(z)  = + +0 *01778-0 .00081  cos  9 +0.00031  cos  20 

+0.00001  cos  39  -0.00001  cos  49  +0.000000  cos  50 
-0.00015  cos  69  +0.03696  sin  9 -0.00353  sin  29 
+0.00641  sin  30  -0.00111  sin  49  +0.00191sin  59-0.00056sin69 
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while  the  actual  duct  coordinates  are  again  obtained  by  integration  of  the  above 
equation,  which  gives  the  result  that 


R-Rl 

BRi 


= 0.09760  loge  (1+0.25^  + 0.10666^-0.00077  sin  {—) 

+ 0.00015  sin  (- j r^L)  -0.00000  sin  {Try)  -0.00000  sin  (^2l) 

+0.00000  sin  -0.00002  sin  (2#>)  +0.05529  [l-cos 

l-cos(£|^ ]+0  .00204  [l-c oe  ( -0 .00026  [l-cos 

(27ry^ 


-0.00168 

+0.00036 


1-cos 


m 


-0.00009  1-cos 


(14) 


If,  on  the  other  hand,  the  formula  for  the  two-dimensional  solution  based 
on  use  of  Ackeret's  formula  is  to  be  sought,  one  has  directly  that 

2.23848  Ai  +3.04608  A 2 = 0.75 
3.66928  A i +3.32288  A2  = 1*125 

provided  it  is  assumed  that  B = 1,  and  0 = 0.25,  and  thus  it  develops  that 
In  the  Normalizing  Mach  No.  Case: 

The  specific  A -values  turn  out  to  be 

Ax  = 0.25  and  A 2 = 0.0625 
and  thus  the  pertinent  slope  function  is  just 

^ W = 1+0  .25 £ + °*0625  1+OJ25I  (13') 

In  the  Small  Conical  Diffusor  Angle  Case 
The  equations  defining  the  A ' s taice  the  evaluated  form 
2.23848  Ai  + 3.04608  A 2 = 0.45 
3.66728  Ai  + 3.32288  A 2 = °*675 


provided  it  is  assumed  that  B = 1.66667  and  0 = 0.15,  and  thus  it  is  found 
that  the  specific  A -values  applicable  to  this  case  are 

Ax  = 0.15  and  A2  = O.C375 

and,  consequently,  the  slope  function  applying  to  the  case  is 

-0*15_  + 0.0375 


(z)  = 


1+0. 25 z 


1+0 .25z 


(14 ») 
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The  actual  radial  coordinates  for  the  duct  giving  minimum  drag  in  this  case 
of  no  increase  in  area  over  that  occupied  by  the  frustum-of-a-cone  basic  shape, 
and  under  the  specific  conditions  of  either  a normalizing  Mach  No.  or  a small 
conical  diffusor  angle,  have  been  calculated  by  means  of  Eqs.  (13)  and  (14)  and 
with  the  integrals  of  Eqs. (13')  and  (14' ),  and  the  final  values  obtained  are 
listed  here  in  Table  XIII, while  the  results  have  also  been  plotted  in  Figs. 7 and  8. 


TABLE  XIII 

DUCT  COORDINATES  FOR  CASE  OF  NO  INCREASE  IN  AREA  OVER  THAT  OCCUPIED 
BY  THE  FRUSTUM-OF-A-CONE  BASIC  SHAPE 


Normalizing  Mach  No.  Case 
B = 1.0  and  p = 0.25 

Small  Conical  Diffusor  Angle  Case 
B = 1.66667  and  p =0.15 

y 

R-R! 

BRi 

R-Rl 

BRi 

By  Use  of  Formula 
(13) 

By  Use  of  Formula 
(13») 

By  Use  of  Formula 
(14) 

By  Use  of  Formula 
(14') 

0 

0 

0 

0 

0 

0.4 

0.09289 

0.1 

0.05575 

0.06 

0.8 

0.19364 

0.2 

0.11619 

0.12 

1.2 

0.29619 

0.3  ' 

0.17773 

0.18 

1.6 

0.40021 

0.4 

0.24013 

0.24 

1.8 

0.45246 

0.45 

0.27147 

0.27 

2.2 

0.55620 

0.55 

0.33373 

0.33 

2.4 

0.60785 

0.6 

0.36470 

0.36 

2.8 

0.70669 

0.7 

0.42398 

0.42 

3.0 

0.75000 

0.75 

0.45000 

0.45 
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duct  contour  for  least  drag  in  case  of  no  increase  in 

AREA  OVER  THAT  OCCUPIED  BY  THE  FRUSTUM -OF ~ A-CONE 
BASIC  SHAPE 

NORMALIZING  MACH  NO  CASE 


FIGURE  7 

DUCT  CONTOUR  FOR  LEAST  DRAG  IN  CASE  OF  NO  INCREASE  IN 
AREA  OVER  THAT  OCCUPIED  BY  THE  FRUSTUM- OF- A- CONE  BASIC 
SHAPE  — SMALL  CONICAL  D/FFUSOR  ANGLE  CASE 


LEGEND < 

MERIDIONAL  CONTOUR  OBTAINED  BY  ASSUMPTION  THAT 

PRESSURE  COEFFICIENTS  ^^F  ARE  RELATED  TO  DUCT 
GEOMETRY  BY  MEANS  OF  EG  (47)  OF  PARTI 

MERIOIONAL  CONTOUR  OBTAINED  BY  ASSUMPTION  THAT 

PRESSURE  COEFFICIENTS  ARE  RELATED  TO  DUCT 

PLANAR  ELEMENTS  THROUGH'  AC  HE  RET’ S TWO-DIMENSION- 
AL FORMULA 
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In  the  small  conical  diffusor  angle  case,  on  the  other  hand,  the  sought 
values  of  A]_  and  A 2 are  determined  by  the  relations 

5.11405  X1  + 4.37629  A 2 = 0.45 

4.96539  A i + 5.06409  A 2 = 0.9756 

so  that  when  the  solution  of  this  pair  of  equations  is  obtained  it  is  seen  that 

Ai  = 0.33399  and  A 2 = 0.13485 

in  this  instance. 

It  follows,  consequently,  that  the  desired  slope  function  is 

1 <■>  ■ %oi!  -0-15495  isfe; +°-01564 

+0.00296  cos  9 +0.00036  cos  29  -0.00007  cos  30 

-0.00001  cos  49+0.00001  cos  5q  -0.00010  cos  69 


+0.03253  sin  9 +0.01303  sin  29  +0.00536  sin  30 
+0.00454  sin  40+0.00162  sin  50  +0.00249  sin  60. 


The  equations  defining  the  actual  external  contour  of  the  duct  giving 
minimum  drag  under  the  above-stated  conditions  then  may  be  written  down,  by 
integration  of  the  slope  functions  just  derived.  When  this  is  done  it  turns 
out  that 

In  the  Normalising  Mach  No.  Case: 


R"R1 


gjjj  = 3.55836  logs  (1+0.25^-0.45258^ 

+0.00231  sin  +0.00027  (^p)  -0.00002  sin 


+0.00091 


+0.00509 

+0.00055 

[l-cos(i^ +0.00034 |l-cos  (2^ 


15) 


while 
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In  the  Small  Conical  Diffusor  Angle  Case: 


R-R 


— i = 5.49524  logg  (l+0.25y)  -0.52568 T 

+0.00285  sin  ( - jf)  +0.00017  sin^^*) -0 .00002  sin  (a'jr) 
+0.00000  sin  (if»)  -0.00000  sin -0.00002  sin  (2 «}) 


3 

+0.03106  1-cos  .00622  [l-cos(^pJ  +0.00171  [l-cos  (^)| 

+0.00108  (l-cos(^p^K) .000 31  [Il-cos^^J  +0.00040  jl-cos  (2»Y)| 


(16) 


On  the  other  hand,  if  the  formula  for  the  two-dimensional  solution,  that 
is  based  on  use  of  Ackeret's  formula  for  relating  the  local  pressures  to  the 
local  duct  slopes,  is  employed,  then  one  finds  that 
In  the  Normalizing  Mach  No.  Case: 

2.23848  Al  +3.04608  A 2 - 0.75 
3.66928  Al  +3.32288  A 2 = 1*4256 
from  whence  it  follows  that 

A-l  = 0.49489  and  A2  = -0*11746 
in  this  instance. 

The  related  slope  function  for  this  example  then  takes  the  form 

z 


n/._»  _ 0.49489 

'P*  - 1+0 .25 i -0'11746 


1+0. 25 z 


and  the  corresponding  actual  radial  coordinates  are  found  by  integrating  the 
slope  function,  to  wit 

= 3.85892  loge  (1+0 .25z)  -0.46984  z.  (15') 

BRi 

In  the  Small  Conical  Diffusor  Angle  Case: 

2.23848  A 1 + 3.04608  A 2 = 0.45 
3.66928  A 1 + 3.32288  A 2 = 0.9756 
and  consequently  the  A values  valid  in  these  circumstances  are 
Ai  = 0.39490  and  /1 2 = -0*14247. 
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The  sought  slope  function  say  then  be  written  down  straight  off  as 

f)  / 2)  — 0.59490  q 14747  z - 
I { 1 “ 1+0 .25 z 1+0.252 

so  that  the  corresponding  definition  of  the  actual  radial  coordinates  comes 
out  to  be 


R— 'R 

= 3.85912  loge  (1+0.25^  - 0.56988*. 


(16') 


The  radial  coordinates  which  apply  to  the  duct  in  question,  which  has 
minimum  drag  in  the  case  where  the  area  enclosed  between  the  sought  duct  contour 
and  the  chord  line  joining  the  and  N2  end-points  is  equal  to  the  area  en- 
closed between  this  Ni  and  N2  chord  line  and  the  parabolic  curve  defined  by 
the  equation 


R _ 
Rl  " 


Y Y X"*  ^1 

1+  0 l + °‘2  £ "°*2  ?T2  * “l 


Rl  ' Ri  Ri“ 

have  been  calculated  by  use  of  Eqs.  (15),  (16),  (15'),  and  (16')  under  the 
specific  conditions  that  a normalizing  Mach  No.  is  being  used  or  that  a small 
conical  diffusor  angle  is  being  treated,  and  the  results  are  entered  into 
Table  XIV  below.  The  results  have  also  been  illustrated  by  means  of  the  plots 
of  Figs.  9 and  10. 
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TABLE  XIV 

DUCT  COORDINATES  FOR  CASE  WHERE  SOUGHT  BEST  SHAPE 
ENCLOSES  SAME  AREA  AS  THAT  COMPRISED  BETWEEN  »1M2 
CHORDLINE  AND  AN  ARBITRARILY  CHOSEN  PARABOLIC  CONTOUR 


w 

Normalizing  Mach  No.  Case 
B=1.0  and  0 =0.25 

Small  Conical  Diffusor  Angle  Case 
B = 1.66667  and  0 = 0.15 

R-Rx 

R-Rx 

BRx 

BRi 

y 

By  Use  of  Formula 

By  Use  of  Formula 

By  Use  of  Formula 

By  Use  of  Formula 

(15) 

(15') 

(16) 

(16') 

0 

o 

0 

0 

0 

0.4 

0.17018 

0.17986 

0.13302 

0.13986 

0.8 

0.32078 

0.32769 

0.24334 

0.24769 

1.2 

0.44614 

0.44862 

0.32768 

0.32862 

1.6 

0.54943 

0.54667 

0.38937 

0.38667 

1.8 

0.59279 

0.58811 

0.41183 

0.40811 

2.2 

0.66482 

0.65752 

0.44233 

0.43752 

2.4 

0.69395 

0.68608 

0.45081 

0.44607 

2.8 

0.73666 

0.73211 

0.45401 

0.45210 

3.0 

0.75000 

0.75000 

0.45000 

0.45000 
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DUCT  CONTOUR  FOR  LEAST  DRAG  IN  CASE  WHERE  THIS  BEST 
DUCT  SHAPE  ENCLOSES  THE  SAME  AREA  AS  THAT  COMPRISED 
BETWEEN  THE  N,  N2  CHORDLINE  AND  THE  PARABOLIC  CURVE 
DEFINED  ACCORDING  TO  THE  EQUATION 

i • w, -°*  w ■ ? 

NORMALIZING  MACH  NO.  CASE 


DUCT  CONTOUR  FOR  LEAST  DRAG  IN  CASE  WHERE  THIS  BEST  SHAPE 


ENCLOSES  THE  SAME  AREA  AS  THAT  COMPRISED  BETWEEN  THE  Nx  Nt 
CHORDLINE  AND  THE  PARABOLIC  CURVE  DEFINED  ACCORDING  TO  THE 
EQUATION : 

SMALL  CONICAL  DIFFUSOR  ANGLE  CASE 


FIGURE  10 


LEGEND < 

— meridional  contour  obtained  by  assumption  that 

PRESSURE  COEFFICIENTS  ARE  RELATED  TO  DUCT 

GEOMETRY  BY  MEANS  OF  EO  {47}  OF  PART  I 

MERIDIONAL  CONTOUR  OBTAINED  BY  ASSUMPTION  THAT 

PRESSURE  COEFFICIENTS  ARE  RELATED  TO  DUCT 

PLANAR  ELEMENTS  THROUGH”  AC  BERET'S  TWO- DIMENSION- 
AL FORMULA 
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DETERMINATION  OF  THE  BEST  DUCT  CONTOUR  UNDER  CONDITION  (a) 
TOGETHER  WITH  CONDITION  (d) 


8.  Calculation  of  the  Coefficients  Hn  and  Kk  for  Perimetral  Condition  on  Volume 
In  the  previous  cases  the  B and  C values  which  were  needed  for  insertion 
into  Eos.  (65)  of  Part  I,  were  evaluated  by  means  of  Eqs.  (64),  but  now  in  this 
instance  where  it  is  required  to  invoke  the  perimetral  condition  having  to  do 
with  the  volume  enclosed  by  the  duct,  as  stated  formally  in  the  (d)  condition  of 
the  Introduction,  it  will  be  necessary  to  make  use  of  the  following  formulations, 
instead  of  the  ones  listed  as  Eos.  (64)  of  Part  I (the  new  coefficients  will  be 
distinguished  from  the  former  ones  by  affixing  an  asterisk  to  the  symbol): 

B*-  f [a, A 2 i, Ce>s  TO 

J’o 

- -2  Jlt  X2’  ' 2^2  ‘ tor  71  odd 


for  71  even 


z r A,  + (l* A 


Likewise 


c;  = j [X,4z\ 4 *$X I+/36I, s/7!  n&d& 

l a. 


•for  71  -0 


for  7?  odd 


~2i; ^ + 2/b£,*. 


X 


f or  71 


even 


for  7?  s 0 
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Consequently,  by  carrying  out  these  coaputations,  for  n running  up  to  6,  the 
following  explicit  values  for  the  sought  coefficients  are  found: 

B0*  = 2.14159  Ax  +5.89049  A 2 

Bi*  = -2.62605  X 2 

B2*  = 0.17905  } 2 

B3*  = -0  .29178  \ 2 

B4*  = 0.04476  A 2 

B5*  = -0 .10504  A 2 

B6*  - 0.01989  A 2 

and 

Cx*  = 2.00000  Ai  +5.66904  A 2 

C2*  = -2.06250  A 2 

C3*  = 0.66667  Ax  +1*35811  A2 

C4*  = -1.05125  A 2 

C5*  = 0.40000  Ai  +0.82155  X2 

C6*  = -0.63750  A 2 

If  the  same  system  of  notation  is  employed  here  as  was  utilized  in  the 
previous  Articles,  it  will  be  seen  that  the  equations  which  determine  the 

w 

values  of  the  IT  and  H coefficients,  which  are  obtained  by  use  of  the  newly 
determined  B and  C constants,  will  take  the  following  explicit  concrete  form  in 
the  present  instance 

0.87169  H0*  +0.00505  H2*  +0.00188  H4*  +0.00100  H6*  = 5.61601 

+0.02551  Li  +0.00440  L3  +0.00152  L5  +0.11599  L0 

0.00190  H0*  +1.02546  H2*  +0.00131  H4*  +0.00074  He*  = 0.09616 

-0.00843  Lx  +0.00795  L3  +0.00157  L5  +0.11599  L2 

0.00100  H0*  +0.00181  H2*  +0-99434  H4*  +0.00054  He*  = -0.04987 

-0.00169  Lx  -0.00566  Lj  +0.00566  L5  +0.11399  L4 

0.00063  H0*  +0.00125  H2*  +0.00065  H4*  +0.095658  H6*  = -0.03177 

-0.00072  Lx  -0.00147  Lg  -0.00299  L5  +0.11399  L6 
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and  likewise 


**  -K#  ** 

0.87169  H0  +0.00505  H2  +0.00188  H4  +0.00100  Hg  = 6.76705 

+0.02531  L{  +0.00440  L3  +0.00132  L5  +0.11399  LQ 

wJt  UJt 

0.00190  H0  +1*02546  H2  +0.00131  H4  +0.00074  Hg  + 0.01325 

-0.00845  +0.00793  Lg  +0.00157  L5  +0.11399  L2 

0.00100  H0**  +0.00181  H2**  +0.99434  H4**  +0.00054  H6**  = -0.05063 
-0.00169  L{  -0.00566  Lg  +0.00366  Lg  +0.11599  L4 
0.00065  H0**  +0.00125  H2**  +C .00065  H4**  +0.95638  H6**  = -0.04244 
-0.00072  L{  -0.00147  Lg  -0.00299  Lg  +0.11399  Lg 
while  ip  like  manner,  the  coefficients  with  odd  subscripts  may  be  deterained  from 
0.96562  H3*  +0.00171  Hg*  +0.00060  H5*  = +0.00901  L2  +0.00251  L4 

+0.00132  Lg  +0.11399  Lx 

0.00295  Hx*  +1.00167  Hg*  +0.00045  H5*  = -0.00540  L2  +0.00537 

+0.00171  Lg  +0.11399  Lg 
0.00197  H-l*  +0.00086  Hg*  +0.99258  Hg*  = -0.00129  L2  -0.00418  L4 

+0.00419  Lg  +0.11399  Lg 

and  likewise 

vw  ww  | 

0.96562  Hx  +0.00171  H +0.00060  Hg  = -2.81963  +0.00901  L2 

+0.00251  L4  +0.00132  Lg  +0.11399  Lx 
0.00295  Hx**  +1.00167  H3**  +0.00045  H5**  = -0.25289  -0.00540  L2 

+0.00537  L4  +0.00171  L^  +0.11399  L3 
0.00197  Hx  +0.00086  H3  +0.99258  H5  = -0.06926  -0.00129  L^ 

-0.00418  L4  +0.00419  Lg  +0.11399  Lg 

Thus,  if  the  L and  L'  values  are  guessed  to  be  zero  for  purposes  of  making 
a first  approximation,  it  will  be  found  that  the  first-stage  solutions  to  these 
equations  will  produce  the  following  results  for  H: 
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H/W  = 4.14828 
H2*^  = -0 .10146 
H/(l)  = -0.05414 


and 


Hi 


*(l) 


h5 


*<*>  * HB#(l)  = 0 


H, 


*(1) 

, - -0.03578 


and  by  use  of  these  values  It  can  be  determined  that  the  corresponding  K values 

» 

are  as  follows 


*(l) 

Kx  * 2.68804 


and 


K 


*(D  _ . *U>  m r *tt>  _ 

j = K4  - Kg  - 0 


Kj*^  = 0.80677 
k5*^  = 0*45582 

and  consequently  the  values  of  L and  L'  may  actually  be  deterained  in  first 
approximation  by  use  of  Eqs.  (66')  of  Part  I,  instead  of  assuming  them  to  be 


zero  as  was  originally  done, 
factors  thus  turn  out  to  be 
L0*  = 1.69114 
L-j*  = -1.39323 
L2*  = -0.34441 
L3*  = 0.04695 
L4*  = -0.17742 
Lg*  = 0.08992 
Lg*  = -0.11272 


The  evaluations  for  the  L*  and  L'*  adjusting 
1* 

and  Lx  * 0.74822 

L^*  s -1-03263 
Lj*  = 0.25727 
L4*  = -0.38338 
Lg*  = 0.15455 
Lg*  = -0.19400 


By  use  of  these  values  of  L#  and  L1#,  a better  second  approximation  for 


the  H's  and  K's  can  then  be  obtained,  and  upon  carrying  out  the  arithmetic 
anew,  it  is  found  that 


*(2) 

*(2) 

Ho 

= 4.39371 

and 

H]/  ' 

* -0.17557 

*(2) 

*(2) 

H2 

= -0.14398 

Hj 

= -0.00904 

a:'2> 

= -0.07678 

H5<2) 

= 0.01280 

«e(2) 

= -0.05075 
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while 


*(2)  *(2) 

Ki  = 2.80785  and  Kg  = 0.47491 

K/U)  = -0.15688  K = -0.02887 

K3*(2)  = 0.84154 

*(2) 

K4  = -0.05647 


In  an  entirely  analogous  way,  the  first  approximation  values  for  the  part 
of  the  H and  K coefficients  that  are  linked  to  the  A.  2 parameter  can  be  obtained 
by  solution  of  the  equations  set  down  above,  and  they  yield  the  following  specific 
values  at  this  stage 


*# 

(1) 

Ho 

= 7, 

.76514 

** 

(1) 

h2 

= -0 

.00146 

(1) 

H4 

= -0 

.05872 

#* 

(1) 

h6 

= -0 

.04945 

Hi 

h5 

**(1) 

H5 


-2.92002 

-0.24587 

-0.06577 


while  the  corresponding  K-values  are  found  to  be 
**(1) 

*1 


**(1) 

Recourse  can  now  be  had  to  the  same  Eqs.  (66')  of  Part  I,  as  employed  pre- 
viously, in  order  to  determine  the  L*#  and  L'#*  factor?  applying  to  this  stage 
in  the  computations;  the  result  of  such  a computation  gives 


**(l) 

= 4.95048 

and 

k4 

- -1.28890 

= -2.70805 

K5**(l) 

= 0.95674 

= 1.64579 

K 

= -0.84177 

it 

.-3 

5.85467 

. T '** 
and 

= 

1.86165 

T **  _ 

L1  * 

-5.55125 

T 1** 

= 

-2.06417 

II 

*C4 

-0.14754 

, 1** 
l5 

* 

0.96828 

** 

L5  = 

0.22216 

T t** 

l4 

= 

-0.78050 

*# 

l4  = 

-0.59442 

- I** 

L5 

s 

0.52456 

1*5  = 

0.50960 

, <** 
Le 

= 

-0.40067 

XX 

L6  = -0.52520 
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and  thus  upon  use  of  these  adjusting  factors  the  second  stage  approximation  for 
the  H**  and  K**  values  may  be  obtained  through  means  of  the  specific  formulae 
set  down  above,  and  it  follows  that  the  new  values,  to  better  approximation,  are 


now 


H, 


H 


H. 


*#(2) 

**(2) 

l 

**(2)  _ 


He 


**(2) 


8 .32478 
-0.02571 
-0.11115 
-0.09304 


and 


(2) 

H1 

= -3.33468 

** 

(2) 

H5 

= -0.21110 

*# 

(2) 

H5 

= -0.02315 

and  likewise  the  second  step  in  the  approximation  to  the  K-coefficients  produces 


the  following  recomputed  values: 
**(2) 

Kx  = 5.22260 
**(2) 

K2  = -3.03105 

^**(2)  _ ^77935 

„ **(2) 

K4  = -1.40629 
**(2) 

K5  = 1.00331 

K6*^(2)  - _o. 90232 


It  is  considered  to  be  sufficiently  accurate  to  stop  the  computations  at 
this  second  stage  in  the  method  of  successive  approximations  because  it  is  only 
the  smallest  (in  absolute  value)  values  of  H and  K which  have  shown  the  greatest 


relative  change  in  going  from  the  first-step  results  to  the  second-step  approxi- 


mations, and  thus  the  ultimate  influence  of  any  such  small  residual  error  upon 
the  final  formula  for  $ ) id.ll  be  scarcely  felt.  If  one  then  divides  the 

sought  slope  function  into  two  parts,  as  was  also  done  previously,  so  that 

')(*)=  fify)  ■<*,  ■'•'UW/U 

then  it  follows  that  the  individual  parts  take  on  the  evaluated  forms  given 
below  when  the  appropriate  values  of  H*,  K*,  H#*,  and  K**,  just  elicited,  have 
been  inserted  into  the  general  formula  for  fjfj)  , which  is  the  counterpart  of 
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Eq.  (62)  of  Part  I that  applied  in  the  previous  cases-  In  the  present  circum- 
stances, therefore,  it  will  be  seen  that 

W =I^255'+  0-12096 

-0-00191  cos  9 +0-00218  cos  29  +0-00003  cos  39 
-0o00009  cos  49  +0  =.00003  cos  59  -0»00100  cos  69 
+0-25142  sin  9 -0-00824  sin  29  +0-04333  sin  39 
-0-00214  sin  49  +0-01292  sin  59  -0-00091  sin  69 

while 

V* = zXi^Sf';  = °-22918 

-0-03635  cos  9 +0-00039  cos  29  +0-00066  cos  39 

-0-00013  cos  49-0,00006  cos  59  -0,00184  cos  69 

+0-46764  sin  9 -0-15926  sin  29  +0-09161  sin  39 

-0-05324  sin  49+0-02729  sin  59  -0,02836  sin  69 

9,  Determination  of  the  /\  l and  A 2 Parameters  for  Perimetral  Condition  on  Volume 
The  stage  has  now  once  again  been  reached  in  the  development  of  the  compu- 
tations where  it  is  possible  to  apply  the  final  specific  perimetral  condition 
which  has  to  do  with  the  restriction  on  the  volume  enclosed  by  the  duct  in  this 
instance-  The  formal  statement  of  this  condition  is  repeated  here  ass 
Case  (d)  In  addition  to  the  closure  condition,  it  is  further  stipulated  now  that 
the  volume  enclosed  between  the  solid  of  revolution  swept  out  by  the  sought  con- 
tour, when  used  as  a generatrix  and  revolved  about  the  duct  axis,  and  the  frustum 
of  a cone  which  is  marked  out  when  the  chordline  Ni  and  N2  is  used  as  a genera- 
trix and  revolved  around  the  same  axis  shall  be  the  equal  to  the  volume  enclosed 
between  the  same  inner  truncated-cone  type  of  body  of  revolution  and  the  one 
which  is  swept  out  as  a result  of  a complete  rotation,  about  the  duct’s  central 
axis,  of  the  parabolic  surface  defined  according  to  the  equation 
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(A) 


in  the  case  of  the  Normalizing  Mach  No..  Mw  = ^2,  or  B = 1,  and  0 = 0.25, 
while  the  generating  curve  is  similarly  taken  to  be 

b f,  * & 16053  bb'b) 

in  the  case  of  the  Small  Conical  Diffusor  Angle,  when  = 1.9456,  B = 1.6667, 
and  0 = 0.15. 


In  addition,  the  parabolic  reference  meridional  line  defined  in  para- 
graph (c)  of  the  Introduction  has  also  been  used  in  applying  tki3  volume-type 
perimetral  condition  for  both  the  Normalizing  Mach  No.  Case  and  the  Small 
Conical  Diffusor  Angle  Case. 

In  all  these  instances  the  first  next  step  is  to  determine  the  appropriate 
values  for  Ai  and  A 2-  This  will  be  done  for  each  of  the  four  "restricted 
volume"  examples  listed  above,  and  they  will  be  treated  in  the  order  set  down  above. 

In  the  Normalizing  Mach  No.  Case, with  (A)  Type  Comparative  Surface 


It  is  found  that  the  specific  form  of  the  comparative  surface  being  used 
in  this  instance  will  produce  the  following  value  for  the  D constant  which 


appears  in  Eq.  (50)  of  Part  Is 

Lz 


D- 


ll 

Wr?  _ 


0.07725 


/-<9*65 


The  corresponding  A -values  may  then  be  determined  from  the  relations 


3.11405  Ai  +5.42245  A2  = 0.75 
7.42077  Ai  + 9.35407  A2  = 1.8468 
which,  upon  solution,  yield  the  values 

A = 0.26992  and  A 2 = -0.01670  . 

The  slope  function  applying  to  this  example  may  then  be  set  down  as 

I)  (.)  = 0.26992  -0.01670  <1+  +0.02882 

1+0.258- 
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-0.00009  cos  9 +0.00058  cos  29  -0.00002  cos  59 
-0.00002  cos  49  +0.00001  cos  59  -0.00024  cos  69 
+0.06005  sin  9 +0.00044  sin  29  +0.01017  sin  39 
+0.00031  sin  49  +0.00303  sin  59  +0.00023  sin  69 
and,  consequently,  the  expression  for  the  actual  radial  coordinates  comes  out 
to  be 


■LSk  = 1.21328  loge  (1+0.25?)  -0.00458?  -0.00417  yz 
BR]_ 


+0.00009  sin  +0.00028  sin  -0.00001  sin  {Try) 

-0.00000  sin^~^)  +0.00000  sin  -0.00004  (2  Vy) 

"l-cos  -0.00021  [l— cos  +o.oo324  [l-cos  (fry) 

(TO0 


+0.05734 

+0.00007 


1-cos 


.00058  l-cos 


(TO 


1+0,00004  [l-cos  (2^ 
(17) 

On  the  other  hand,  one  may  proceed  on  in  a like  manner,  except  for  a 
slight  change  in  the  value  of  D,  in  order  to  obtain  the  ordinates  for  the 
Small  Conical  Diffuser  Angle  Case,  with  (B)  Type  Comparative  Surface 

Here  the  determinative  equation  giving  the  sought  value  of  D is 

[f  . ««7  ■ 

. ) 

The  corresponding  values  for  the  /j  s are  found,  therefore,  from  the 

equations  . 

3.11405  A1  +5,42245  A 2 = °“45 

7.42077  Ai  +9.35407  2 = 1.3437 

so  that  the  specific  values  of  A applying  in  this  example  are 

and  A 2 = -0.07606  . 


X=  0.27695 

Thus  one  may  write  down  the  appropriate  slope  function  almost  immediately 

fll 25£) 


as  (J  (?)  = 0.27695 


1+0.25S 


-0.07606z' 


Ll 


1+0 .25r 


+0.01607 
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+0.00224  cos  9 +0.00057  cos  20  -0.00006  cos  39 
-0.00002  cos  49  +0.00001  cos  59  -0.00014  cos  60 
+0.03406  sin  0 +0.00983  sin  20  +0.00503  sin  39 
+0.00346  sin  49  +0.00150  sin  50  +0.00190  sin  60 


and  it  follows,  upon  integration,  that  the  radial  coordinates  of  the  duct  contour 
in  this  specific  case  are  given  by  the  expression 


R-Rl 

BRi 


1.71628  loge  (1+0.25^  -0.13605y  -0.01901  y2 
+0.00214  sin  (-^J+0 .00027  sin ) -0 .00002  sin  (r  fr) 
-0.00000  sin (^p~)+0 .00000  sin  ^-^i-J-0. 00002  sin  (2JTJ  ) 
+0.03252  l-cos^-^J+O .00469  [l-cos  (^p^jJ+0 .00160  1-cos 
+0.00083  1-cos  ^pJj+0. 00029  |l-cos(^1-J]+0. 00030  1-cos 


(2  *r\ 


(18) 


It  is  again  instructive  to  compare  these  results  with  the  contour  which 
would  result  if  the  two-dimensional  approach  were  to  be  followed, which  relies  on 
use  of  Ackeret's  formula  for  relating  the  local  pressures  to  the  local  duct  slopes. 
Such  determinations  are  easily  carried  out,  following  the  procedure  amply  illus- 
trated previously,  so  that  it  is  at  once  found  that,  for 
The  Normalizing  Mach  No.  Case 

The  determinative  equations  for  the  A <+>  are  this  time 

2.23848  Al  + 3.77304  A2  = 0.75 

5.46072  Xl  + 6.09378  A 2 = 1.8468 

from  whence  it  follows  that 

Al  = 0.34435  and  A 2 = -0.00552  . 

The  corresponding  slope  function  for  the  "two-dimensional"  case  then 

turns  out  to  be  / 0 25  A 

n *(l+  -7—  ) 

r)  (*)  = 0.34435  jj0;2S-z.  -0.00552  1+0.25*- 

and  thus  the  sought  expression  for  the  radial  coordinates  of  the  duct  is 
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Likewise 


R-Rl 

BRi 


= 1 o42156  loge  (l+0.25z)  -0 .01104  * -0.00138^ 


(17 ») 


In  the  Small  Conical  Diffusor  Angle  Case 

The  appropriate  detenainative  equations  for  the  A J are  now 
2.25848  Ai  +3.77504  /l 2 = 0.45 
5.46072  Al  +6.09378  A 2 = 1.3437 

so  that  thus 

Ai  = 0.33431  and  A 2 = -0.07907. 

The  appropriate  slope  function  is  obtained  by  substitution  of  these  values  Ai 
A 2 Into  the  usual  general  expression,  with  the  result  this  time  that  one  obtains 
the  elope  as 

<z)  = 0.55451  -°‘ 07907  “ 

with  the  consequence  that  the  "two-dimensional"  duct  of  best  shape  comes  out  to 
be  in  this  case 

R-Rl  o 

= 1.96980  loge  (1+0 .25 7) -0.15814 7 -0.01977  tT . 

The  radial  coordinates  which  apply  to  the  ducts  in  quest! on, which  have 
a minimum  drag  in  the  cases  where  the  volume  enclosed  between  the  solid  of 
revolution  swept  out  by  the  sought  contour  when  revolved  about  the  duct's  axis 
and  the  frustum  of  a cone  which  has  the  chordline  Ni  N2  as  a generatrix  is 
equal  to  the  volume  enclosed  between  this  same  inner  truncated-cone  type  of 
body  of  revolution  and  the  one  which  is  swept  out  as  a result  of  a complete 
rotation,  about  the  duct's  central  axis,  of  the  parabolic  curve  defined  accord- 
ing to  either  of  the  two  equations  given  in  the  (d)  part  of  the  Introduction, 
have  been  calculated  by  use  of  Eqs.  (17),  (18),  (17 ' ) , and  (18')  under  the 
specific  condition  that  a normalizing  Mach  No.  is  being  used  in  one  case, 
while  that  a small  conical  diffusor  angle  is  being  assumed  in  the  other,  and 
the  results  are  entered  into  Table  XV  below.  The  results  have  also  been 
illustrated  by  means  of  the  plots  of  Figs.  (11)  and  (12). 


(1+  2^ ,) 

1+0 .25  z 
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DUCT  CONTOUR  FOR  LEAST  DRAG  IN  CASE  WHERE  TH/S  BEST 
SHAPE  ENCLOSES  THE  SAME  VOLUME  AS  THAT  COMPRISED 
BETWEEN  THE  INNER  FRUSTUM  - OF  - A-  CONE  SURFACE  SWEPT 
OUT  BY  THE  N,Nt  CHORDLINE  AND  THE  OUTER  SURFACE  SWEPT 
OUT  BY  A PARABOLIC  CURVE  HAVING  THE  EQUATION 

%•'**%*  007725  (ww%) 

NORM  A LIZ  (NS  MACH  HO.  CASE 


DUCT  CONTOUR  FOR  LEAST  DRAG  IN  CASE  WHERE  THIS  BEST 
SHAPE  ENCLOSES  THE  SAME  VOLUME  AS  THAT  COMPRISED  BE- 
TWEEN THE  INNER  FRUSTUM  - OF- A- CONE  SURFACE  SWEPT  OUT  BY 
THE  Nj  N2  CHORDUNE  AND  THE  OUTER  SURFACE  SWEPT  OUT  BY  A 
PARABOLIC  CURVE  HAVING  THE  EQUATION * 

' ***]%  *OJ6038  (•%/••£  -ft  ) 

- SMALL  CONC/AL  DIFFUSOR  ANCLE  CASE  * 


meridional  contour  obtained  by  assumption  that 

PRESSURE  COEFFICIENTS  ARE  RELATED  TO  DUCT 

GEOMETRY  BY  MEANS  OF  EO  (47)  OF  PART  1 
’ MERIDIONAL  CONTOUR  OBTAINED  BY  ASSUMPTION  THAT 
PRESSURE  COEFFICIENTS  gfc  ARE  RELATED  TO  DUCT 
PLANAR  ELEMENTS  THROUGH  ACKERET'S  TWO-DIMENSION- 
AL FORMULA 
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TABLE  XV 


DUCT  COORDINATES  FOR  CASE  WHERE  SOUGHT  BEST  DUCT  SHAPE  ENCLOSES 
SAME  VOLUME  AS  THAT  COMPRISED  BETWEEN  INNER  FRUSTUM  OF  A CONE 
SURFACE  SVEPTED  OUT  BY  THE  $1*2  CHORDLINE  AND  THE  OUTER  SURFACE 
SWEPT'  OUT  BY  THE  ARBITRARILY  CHOSEN  PARABOLIC  CONTOURS  DEFINED 
IN  CASE  (d)  OF  INTRODUCTION 


Normalizing  Mach  No.  Case 

Small  Conical  Diffusor  Angle  Case 

p = 1.0  and  p = 0.25 

P = 1.66667  and  P^=  0.15 

R-R]_ 

R-Rx 

y 

BRX 

BRl 

By  Use  of  Formula 

By  Use  of  Formula 

By  Use  of  Formula 

By  Use  of  Formula 

(17) 

(17') 

(18) 

(18') 

0 

0 

0 

0 

0 

0.4 

0.12163 

0.13085 

0.11452 

0.12132 

0.8 

0.24135 

0.24946 

0.21489 

0.21997 

1.2 

0.35297 

0.35772 

0.29647 

0.29856 

1.6 

0.45723 

0.45711 

0.36059 

0.35914 

1.8 

0.50631 

0.50385 

0.38585 

0.38319 

2.2 

0.59822 

0.59203 

0.42389 

0.41967 

2.4 

0.64130 

0.63369 

0.43682 

0.43239 

2.8 

0.71845 

0.71259 

0.44962 

0.44745 

3.0 

0.75000 

0.75000 

0.45000 

0.45000 

Since  it  may  possibly  be  of  more  advantage,  for  comparative  purposes,  to 
to  use  the  same  shape  for  the  "control"  surface,  regardless  of  whether  the  restric- 
tion of  the  perimetral  condi tion  has  to  do  with  the  area  enclosed  in  one  meri- 
dional plane  or  the  volume  enclosed  by  the  revolved  curve,  the  "volume"  condition 
will  also  now  be  applied  to  the  case  where  the  comparative  contour  is  given  by 
the  parabolic  equation  set  down  in  the  (c)  part  of  the  Intro ductlon.  Carrying 
through  the  same  processes,  as  just  illustrated  above,  for  this  new  case,  there- 
fore, it  will  be  necessary  first  of  all  to  find  a new  value  for  D.  This  will  be 
done  for  the  usual  specific  instances  of  the  normalizing  Mach  No-  case  and  the 
small  conical  diffusor  angle  case: 
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D-  B 'R 7 [if  + ^ Z ^/00°  P°r  B^l.0  m4  J*aZS 

*0=  gj^- [f+  ^£-+($-+12 : /'42S<)  *r  B‘,(Aiii ind  Aa/5- 

Consequently  the  rest  of  the  computation  goes  through  in  the  same  straight- 
forward manner  employed  previously,  so  that 
In  the  Normalizing  Mach  No.  Case 

The  appropriate  values  of  A will  be  found  from  solution  of  the  equations 


3.11405  \1  +5.42245  A2  = 0.7500 
7.42077  Ai  +9.35407  X2  = 2.1000 
so  that  the  sought  values  may  be  seen  to  be 

X i = 0.39351  . and  A2  = -0.08768  . 

By  use  of  these  values  of  A in  the  general  expression  for  the  slope  function, 
therefore,  it  will  be  seen  that  in  this  case 

0-25  \ 

n (?)  = 0.39351  ——-0.08768  2^ — Z +0.02750 

1 1+0 .252  1+0. 25 z 

+0.00244  cos  0 +0.00082  cos  20  -0.00007  cos  30 
-0 .00003  cos  40+0.00002  cos  50  -0.00023  cos  60 


+0.05793  sin  0 +0.01072  sin  20  +0.00902  sin  30 
+0.00383  sin  40+0.00269  sin  50  +0.00213  sin  60 


and  then  the  corresponding  radial  coordinate  function  is  found  to  be 

R — Rt  — 

= 2.27548  loge  (1+0.25$  -0.14786^-0.02192^ 

+0.00233  sin  (-^9+0.00039  sin -0-00002  sin  {Try) 

-0.00001  sin^l)+0. 00000  sin^j^-J  -0.00004  sin  {2TTy) 

+0  .05532  [l--cos  P£>]  +0.00512  [l-cos(^9)+0  .00287  [l-cos  (r  y)| 

+0 .00091  [l-cos  Ji|3()jf0  .00051  [l-cos^-^JJ+O .00034  [l-cos  (2^y )| 


(19) 
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Ip  the  Small  Conical  Diffnsor  Angle  Case 

The  appropriate  values  of  A will  be  obtained  from  solution  of  the 
simultaneous  set 

5.11405  A l + 5.42245  A2  = 0.4500 
7.42077  A 1 + 9.35407  A2  = 1.4250 

and  consequently 

Ax  = 0.31663  and  A 2 = “0*09885 

with  the  result  that  the  sought  slope  function  will  be  given  in  the  form 

n(*)  = 0.51663  * — ^ -0*09885  • ICkt 7^1.  1+q .01564 
I l+o.^&sr  1+0 .25 z 

+0.00299  cos  6 +0.00065  cos  20  -0.00007  cos  30 


-0.00002  cos  40+0.00001  cos  50  -0.00013  cos  60 


+0.03338  sin  0+0.01313  sin  20  +0.00466  sin  30 
+0.00458  sin  40+0.00139  sin  50  +0.00251  sin  60  . 

Upon  integration  of  this  function  the  final  equation,  describing  the 
way  in  which  the  external  contour  must  vary  in  this  instance  to  give  least 


drag,  is  obtained,  and  this  computation  may  be  readily  performed  to  give 
H— Ri 

= 2.05732  loge  (1+0 .25 Jr)  -0 .18206 y -0.02471  y 2 

+0.00285  sin  l^~)  +0.00031  sin^^^-J-O .00002  sin  {fy) 

-0.00000  sin (Af*)  +0.00000  sin  (-—^-0 .00002  sin  (2  7T}) 
+0.03188  l-cos(-V J +0.00627  [l-cos^^pJ+O  .00148  [l-cos  {T}  )] 
+0.00109  1— cos(^i‘  ] +0.00026  [l-cos(^J^J  +0.00040 [l-cos  (2 Tfy^ 

(20) 

By  application  of  the  "two-dimensional"  treatment  in  these  cases,  however,  it 


will  be  seen  that 
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In  the  Wo  mail  zing  Mach  No.  Case 

2 .23848  A x 4-5.77504  ^ 2 = 0.75 
5.46072  A i 4-6.09578  A 2 = 2.1006 


and,  hence,  one  gets  that 

Ax  = 0.48189  and  A2  = -0.08712  . 

The  corresponding  slope  function  then  takes  the  for* 

n z(l+  “ 2 *' 

fj  (z)  = 0.48189  i+oT25z  l+0.25z 

and  the  related  value  for  the  actual  radial  coordinates  then  turns  out  to  be 


= 2.62452  loge  (1+0 .25 z)  -0 .17424 z -0 .02178 Z2 


(19») 


Likewise, 

In  the  Small  Conical  Dif fusor  Angle  Case 
One  has  this  time  that 


2.23848  Ax  + 3.77304  A2  = 0.45 
5.46072  A 2 4-  6.09378  A 2 = 1.4256 
and  thus  the  appropriate  values  of  A in  this  case  are 
Ax  = 0.37869  and  A2  = -0.10540  , 

By  substitution  of  these  values  into  the  general  expression  for  the 
slope  function,  therefore,  one  finds  that  now 


0-25 
z(l+~2~ 


z) 


0.37869  1+OJ25Z  -0.10540  1+0j25z 


and  thus  the  sought  "two-dimensional"  solution  for  the  best  duct  in  this  case 
is 


R-R^ 

"brJ 


= 2.35796  log.  ( 1+0 J25 z) -0 .21080 z-0 .02635 z 


(20  •) 
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The  radial  coordinates  which  apply  to  the  duct  which  has  a minimum 
drag  under  the  circumstances  that  the  volume  enclosed  between  the  solid  of 
revolution  swept  out  by  this  best  contour  when  revolved  about  its  axis  and 
the  frustum  of  a cone  which  has  the  chordline  N2N2  as  a generatrix  is  equal 
to  the  volume  enclosed  between  this  same  inner  truncated-cone  type  of  body 
of  revolution  and  the  one  which  is  swept  out  as  a result  of  a complete  rota- 
tion, about  the  duct's  axis,  of  the  parabolic  line  defined  by  means  of  the 
equation  2 

-§r  '^itazf,-az 

in  any  meridional  plane  have  been  computed  by  use  of  Eqs.  (19),  (20), (19'),  and 
(20*)  under  the  specific  condition  that  a normalizing  Mach  No*  is  governing 
the  general  flow  as  well  as  in  the  case  that  a small  conical  diffusor  angle 
type  of  duct  is  under  study.  The  results  have  been  entered  into  Table  XVI 
and  the  results  have  also  been  illustrated  by  means  of  the  plots  constitut- 
ing the  appended  Figs.  (13)  and  (14). 

TABLE  XVI 

DUCT  COORDINATES  FOR  CASE  WHERE  SOUGHT  BEST  DUCT  SHAPE  IS  ONE  WHICH 
ENCLOSES  SAME  VOLUME  AS  THAT  LEFT  BETWEEU  INNER  FRUSTUM  OF  A CONE 
SWEPT  OUT  BY  NlNo  CHORDLINE  AND  THE  OUTER  SURFACE  SWEPT  OUT  BY  THE 
SAME  PARABOLIC  CONTOUR  USED  AS  A REFERENCE  IN  THE  AREA-RESTRICTION  CASE 


Normalizing  Mach  No.  Case 

Small  Conical  Diffusor  Angle  Case 

B = 1.0  p 

= 0.25 

B = 1.6667 

p = 0.15 

R-Fa 

R-Rx 

y 

BRi 

BRX 

By  Use  of  Formula 

By  Use  of  Fbimula 

Ey  Use  of  Formula 

By  Use  of  Formula 

(19) 

(19') 

(20) 

(20') 

0 

0 

0 

0 

0 

0.4 

0.16625 

0.17697 

0.12883 

0.15620 

0.8 

0.31665 

0.52517 

0.25903 

0.24440 

1.2 

0.44397 

0.44812 

0.32567 

0.32775 

1.6 

0.54989 

0.54855 

0.59032 

0.33865 

1.8 

0.59444 

0.59097 

0.41411 

0.41131 

2.2 

0.66797 

0.66145 

0.44625 

0.44208 

2.4 

0.69717 

0.68990 

0.45474 

0.45054 

2.8 

0.75850 

0.73402 

0.45596 

0.45458 

3.0 

0.75000 

0.75000 

0.45000 

0.45000 

art 


DUCT  CONTOUR  FOR  LEAST  DRAG  IN  CASE  WHERE  THIS  BEST 
SHAPE  ENCLOSES  THE  SAME  VOLUME  AS  THAT  LEFT  BETWEEN 


THE  INNER  FRUSTUM -OF- A - CONE  SWEPT  OUT  BY  THE  Nf  N£ 
CHORDLINE  AND  THE  OUTER  SURFACE  SWEPT  OUT  BY  THE  SAME 
PARABOLIC  CONTOUR  USED  AS  A REFERENCE  IN  THE  AREA- 
RESTRICTED  CASE,  AND  THE  EQUATION  OF  WHICH  IS> 

t,'  '*■*'  n; ' °x w; * 02 ** ' t 

NORMALIZING  MACH  NO  CASE 


DUCT  CONTOUR  FOR  LEAST  DRAG  IN  CASE  WHERE  THIS  BEST 
SHAPE  ENCLOSES  THE  SAME  VOLUME  AS  THAT  LEFT  BETWEEN 
THE  INNER  FRUSTUM  -OF- A -CONE  SWEPT  OUT  BY  THE  NfNe 
CHORDLINE  AND  THE  OUTER  SURFACE  SWEPT  OUT  BY  THE 
SAME  PARABOLIC  CONTOUR  USED  AS  A REFERENCE  IN  THE 
AREA-  RESTRICTION  CASE,  AND  THE  EQUATION  OF  WHICH  IS: 


SMALL  CONICAL  D/FFUSOR  ANGLE  CASE 


LEGEND' 

— meridional  contour  OBTAINED  BY  ASSUMPTION  that 
PRESSURE  COEFFICIENTS  f'fc  ARE  RELATED  TO  DUCT 
GEOMETRY  BY  MEANS  Of  ~EQ.  ( 47 ) OF  PART  X 

MERIDIONAL  CONTOUR  OBTAINED  BY  ASSUMPTION  THAT 

PRESSURE  COEFFICIENTS  ARE  RELATED  TO  DUCT 

PLANAR  ELEMENTS  THROUGH  ACKERET'S  TWO-DIMENSION- 
AL FORMULA 
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Conclusions 


10. 


By  reference  to  the  above-given  tabular  results  for  any  of  the  specific 


sets  of  circumstances  considered,  or  by  visual  comparison  in  the  individual 

R-Rl 

figures,  it  may  be  seen  that  in  all  cases  the  values  of 


(the  radial  co- 


ordinates of  the  sought  duct)  when  calculated  through  means  of  the  formula 
which  was  derived  in  Part  I as  Eq»  (47),  for  relating  the  geometrical  features 
of  the  duct  to  the  pressures  supported  by  its  surface,  do  not  differ  in  any 
sizeable  degree,  to  engineering  standards  of  precision,  from  the  corresponding 
duct  coordinates  which  are  obtained  by  mere  application  of  the  "two-dimensional" 
procedure  for  determination  cf  the  optimal  duct  shape,  which  depends  on  the  use 


of  Ackeret's  formula  for  relating  the  local  pressure  coefficient. 


P-P. 


22. 


P U 

eo  CO 


2 > to 


the  local  slopes  of  the  constraining  surface. 

As  has  been  pointed  out  in  the  Introduction,  this  result  should  not  be 
interpreted  as  being  of  inconsecuential  import,  because  it  nevertheless,  de- 
spite the  similarity  in  contours,  is  going  to  be  found  true  that  the  actual 

p-p 

values  of  5 that  are  obtained  by  use  of  Ecu  (47)  of  Fart  I for  relating 

P U jf 

*00  OO 


the  local  planar  elements  of  the  ducts  in  Question  to  the  pressures  will  be 
appreciably  different  from  the  corresponding  values  for  the  pressures  which 
are  obtained  by  use  of  the  "two-dimensional"  treatment.  The  extent  and  import 
of  these  differences  in  pressures,  which  arise  under  the  dictates  of  these  two 
contrasting  assumptions,  may  most  simoly  be  realized  by  reference  to  the  plots 
of  the  longitudinal  components  of  the  induced  velocities,  produced  in  these 
two  cases,  which  are  diagrammed  in  Figs.  3 and  4 for  each  of  the  harmonic  com- 
ponents making  up  the  total  induced  velocity;  it  will  be  seen  that  for  the 
lower  values  of  n and  k there  are  marked  differences  in  the  plotted  curves 
for  the  longitudinal  velocity  components. 
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On  the  other  hand,  if  one  is  merely  interested  in  the  determination  of  the 
ordinates  of  the  best  duct,  without  concern  about  how  to  compute  the  correspond- 
ing pressures,  it  is  obvious  that  this  general  result,  adduced  from  the  compari- 
sons now  carried  out  and  on  the  basis  of  the  assurances  given  in  the  text  about 
the  wide  generality  of  the  illustrative  examples  selected  for  computation,  may 
be  of  quite  wide  practical  interest,  because  the  solution  of  the  problem,  as 
far  as  determining  the  optimal  duct  ordinates  is  concerned,  may  now  be  attained 
by  use  of  the  simplest  of  procedures,  depending  merely  upon  application  of 
Ackeret's  two-dimensional  formula.  Of  course,  it  will  still  be  necessary  to 
use  the  more  complicated  formulae  developed  in  Part  I for  relating  the  duct  con- 
tours to  the  pressures  exerted  on  its  external  walls  if  one  wishes  to  obtain 
suitably  exact  values  for  the  pressures. 
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